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We present the statistical properties of the first version of the Cold Core Catalogue of Planck Objects (C3PO), in terms of their spatial distribution, 
' temperature, distance, mass, and morphology. We also describe the statistics of the Early Cold Core Catalogue (ECC) that is a subset of the 
complete catalogue, and that contains only the 915 most reliable detections. ECC is delivered as a part of the Early Release Compact Source 
Catalogue (ERCSC). We have used the CoCoCoDeT algorithm to extract about 10 thousand cold sources. The method uses the IRAS lOOyura data 
as a warm template that is extrapolated to the Planck bands and subtracted from the signal, leading to a detection of the cold residual emission. 
We have used cross-correlation with ancillary data to increase the reliability of our sample, and to derive other key properties such as distance and 
mass. 

Temperature and dust emission spectral index values are derived using the fluxes in the IRAS 100 fim band and the three highest frequency Planck 
bands. The range of temperatures explored by the catalogue spans from 7 K to 17 K, and peaks around 13 K. Data are not consistent with a constant 
value of the associated spectral index f> over the all temperature range./? ranges from 1.4 to 2.8 with a mean value around 2.1, and several possible 
k> | scenarios are possible, including P(T) and the effect of multiple temperature components folded into the measurements. 

' i . For one third of the objects the distances are obtained using various methods such as the extinction signature, or the association with known 
■ molecular complexes or Infra-Red Dark Clouds. Most of the detections are within 2 kpc in the Solar neighbourhood, but a few are at distances 
greater than 4 kpc. The cores are distributed over the whole range of longitude and latitude, from the deep Galactic plane, despite the confusion, to 
high latitudes (> 30°). The associated mass estimates derived from dust emission range from 1 to 10 5 solar masses. Using their physical properties 
such as temperature, mass, luminosity, density and size, these cold sources are shown to be cold clumps, defined as the intermediate cold sub- 
structures between clouds and cores. These cold clumps are not isolated but mostly organized in filaments associated with molecular clouds. The 
Cold Core Catalogue of Planck Objects (C3PO) is the first unbiased all-sky catalogue of cold compact objects and contains 10783 objects. It gives 
an unprecedented statistical view to the properties of these potential pre-stellar clumps and offers a unique possibility for their classification in 
terms of their intrinsic properties and environment. 

Key words. Cold Cores, Galaxy, Source extraction 
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The main difficulty in understanding star formation lies in the 
vast range of scales involved in the process. If star formation it- 
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self is the outcome of gravitational instability occurring in cold 
and dense structures at sub-parsec scales, the characteristics of 
these structures (usually called pre-stellar cores) depend on their 
large-scale environment, up to Galactic scales because their for- 
mation and evolution is driven by a complex coupling of self- 
gravity with cooling processes, turbulence and magnetic fields, 
to name a few. To progress in the understanding of star forma- 
tion pre-stellar cores need to be observed, in a variety of environ- 
ments. More importantly, broad surveys are required to address 
statistical issues, and probe theoretical predictions regarding the 
initial mass function (IMF) largely determined at the stage of 
fragmentation of pre-stellar cores. 

Unfortunately, the properties of the pre-stellar cores are still 
poorly known mostly because of observational difficulties. The 
total number of Galactic pre-stellar cores is estimated to be 
around 3 x 1 5 (Clemens et al. 1 99 1 ) but most of them have so far 
escaped detection, simply because they are cold and immersed 
in warmer (therefore brighter) environments. 

The thermal dust emission of nearby molecular clouds has 
been mapped from the ground in the millimeter and submillime- 
ter ranges with instruments such as SCUBA, MAMBO, SIMBA, 
and Laboca. Because of limited sensitivity, but also the presence 
of the atmospheric fluctuations that call for beam-throw of at 
most a few arcmin, the studies have concentrated on the bright- 
est and most compact regions that are already in an active phase 
of star formation. Thanks to sub-arcminute resolution, these ob- 
servations (together with dedicated molecular line studies) have 
been the main source of information also on the structure of 
the pre-stellar cores (Motte et al. 1998; Curtis & Richer 2010; 
Hatchell et al. 2005; Enoch et al. 2006; Kauffmann et al. 2008). 

Many compact clouds were detected as absorption features 
on photographic plates. A new population of thousands of cold 
dark clouds was discovered by observations of mid-infrared 
absorption towards the bright Galactic background (MSX and 
ISOGAL surveys; see Egan et al. 1998; Perault et al. 1996). The 
absorption studies are, however, strongly biased towards the low 
latitudes and do not directly provide information on the tem- 
perature of the detected sources. For a definitive study of the 
cold cloud cores, one must turn to high resolution observations 
in the submillimetre or millimetre range (Andre et al. 2000). The 
Bolocam Galactic Plane Survey (BGPS) is producing mm data 
for the central part of the Galactic plane (Aguirre et al. 2010). 
The first results suggest that at kpc distances, even with a half ar- 
cmin resolution, one is detecting mainly cluster forming clumps 
rather than cores that would produce, at most, a small multiple 
system (Dunham et al. 2010). 

Balloon borne experiments have provided larger blind sur- 
veys of higher latitudes. PRONAOS discovered cold condensa- 
tions also in cirrus-type clouds (Bernard et al. 1999; Dupac et al. 
2003) Similarly, Archeops (Desert et al. 2008) detected hundreds 
of sources with temperatures down to 7 K. The latest addition 
to the balloon borne surveys is the BLAST experiment which 
has located several hundred submillimetre sources in Vulpecula 
(Chapin et al. 2008) and Vela (Netterfield et al. 2009; Olmi et al. 

2009) , including a number of cold and probably pre-stellar cores. 
Since its launch in May 2009, the Herschel satellite has al- 
ready provided hundreds of new detections of both starless and 
protostellar cores (Andre et al. 2010; Bontemps et al. 2010; 
Konyves et al. 2010; Molinari et al. 2010; Ward-Thompson et al. 

2010) . There is an intriguing similarity between the core mass 
function (CMF) derived from these data, and the IMF that need 
to be investigated in different environments, towards the inner 
Galaxy in particular. The Herschel studies will eventually cover 
a significant fraction of the Galactic mid-plane and the central 



parts of the nearby star-forming clouds but cannot cover high 
Galactic latitudes where star formation is known to occur. In this 
endeavor the main challenge is how to locate the cores because, 
even with Herschel, detailed studies must be limited to a small 
fraction of the whole sky. 

The Planck 1 satellite (Tauber et al. 2010) improves 
over the previous studies by providing an all-sky submillime- 
tre/millimetre survey that has both the sensitivity and resolution 
needed for the detection of compact sources. The shortest wave- 
length channels of Planck cover the wavelengths around and 
longwards of the intensity maximum of the cold dust emission: 
v 2 B v (T = 10K) peaks close to 300 pm while, with a temperature 
of T ~ 6 K, the coldest dust inside the cores has its maximum 
close to 500 pm. Combined with far-infrared data such as the 
IRAS survey, the data enable accurate determination of both the 
dust temperature and the spectral index. We use the Planck ob- 
servations to search for Galactic cold cores, i.e. compact cloud 
cores with colour temperatures below 14 K. Because of the lim- 
ited resolution, we are likely to detect mainly larger clumps in- 
side which the cores are located. The cores will be pre-stellar 
objects before (or at the very initial stages) of the protostel- 
lar collapse, or possibly more evolved sources that still contain 
significant amounts of cold dust. The Cold Core Catalogue of 
Planck Objects (C3PO) which will be made public at the end 
of the Planck proprietary period, will be the first all-sky cat- 
alogue of cold cloud cores and clumps. It will reveal the lo- 
cations where the next generations of stars will be born and 
will provide an opportunity to address a number of key ques- 
tions related to Galactic star formation: What are the character- 
istics of this source population? How does the distribution of the 
cores/clumps correlate with the current star formation activity 
and the location of the molecular cloud rings and the spiral arms? 
How are the sources related to large-scale structures like the FIR 
loops, bubbles, shells, and filaments? Are there pre-stellar cores 
at high latitudes? How much do the core properties depend on 
their environment? Investigations such as these will help us un- 
derstand the origin of the pre-stellar cores, the instabilities that 
initiate the collapse, and the roles of turbulence and magnetic 
fields. The catalogue will prove invaluable for follow-up studies 
to investigate in detail the internal properties of the individual 
sources. 

In this paper we describe the general properties of the 
current cold cores catalogue that is based on data that the 
Planck satellite has gathered during its first two scans of the 
full sky. In particular, we will describe the statistics of the 
Early Cold Cores Catalogue (ECC) that is part of the recently 
published Planck Early Release Compact Source Catalogue 
(ERCSC Planck Collaboration 2011c). ECC forms a subset of 
the full C3PO and contains only the most secure detections of all 
the sources with colour temperatures below 14 K. The final ver- 
sion of C3PO will be published in 2013. For historical reasons, 
we use "Cold Cores" to designate the entries in the C3PO and in 
the ECC, and similarly in much of this paper. However, as this 
paper and the companion paper (Planck Collaboration 201 lr, 
hereafter Paper II) demonstrate, most of these are more correctly 
described as "cold clumps", intermediate in their structure and 



1 Planck (http://www.esa.inl/Planck) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries 
France and Italy), with contributions from NASA (USA) and telescope 
reflectors provided by a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 
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physical scale between a true pre-stellar core and a molecular 
cloud. 

Planck (Tauber et al. 2010; Planck Collaboration 2011a) is 
the third generation space mission to measure the anisotropy of 
the cosmic microwave background (CMB). It observes the sky 
in nine frequency bands covering 30-857 GHz with high sensi- 
tivity and angular resolution from 31' to 5'. The Low Frequency 
Instrument LFI; (Mandolesi et al. 2010; Bersanelli et al. 2010; 
Mennella et al. 201 1) covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HFI; 
Lamarre et al. 2010; Planck HFI Core Team 2011a) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolometers 
cooled to 0. 1 K. Polarization is measured in all but the highest 
two bands (Leahy et al. 2010; Rosset et al. 2010). A combina- 
tion of radiative cooling and three mechanical coolers produces 
the temperatures needed for the detectors and optics (Planck 
Collaboration 2011b). Two Data Processing Centers (DPCs) 
check and calibrate the data and make maps of the sky (Planck 
HFI Core Team 2011b; Zacchei et al. 2011). Planck's sensitiv- 
ity, angular resolution, and frequency coverage make it a pow- 
erful instrument for galactic and extragalactic astrophysics as 
well as cosmology. Early astrophysics results are given in Planck 
Collaboration, 201 lh-z. 

2. Source Extraction 

2.1. Data Set 

As cold cores are traced by their cold dust emission in the 
submillimetric bands, we use Planck channel maps of the HFI 
at 3 frequencies : 353, 545 and 857 GHz as described in de- 
tail in Planck HFI Core Team (2011b). The temperature maps 
at these frequencies are based on the first two sky surveys 
of Planck, provided in Healpix format (Gorski et al. 2005) at 
nside=2048. We give here a very brief summary of the data re- 
duction, cf Planck HFI Core Team (2011b) for further details. 
Raw data are first processed to produce cleaned timelines (TOI) 
and associated flags identifying various systematic effects. The 
data analysis includes application of a low-pass filter, removal 
and correction of glitches, conversion to absorbed power and 
decorrelation of thermal stage fluctuations. For the cold core de- 
tection, and more generally for source detection, Solar System 
objects (SSO) are identified in the TOI data using the publicly 
available Horizon ephemerides and an SSO flag is created to en- 
sure that they are not projected onto the sky. 

Focal plane reconstruction and beam-shape estimates are ob- 
tained using observations of Mars. Beams are described by an 
elliptical Gaussian parameterisation leading to FWHM 6s, given 
in Table 2 of Planck HFI Core Team (2011b). The attitude of 
the satellite as a function of time is provided by the two star 
trackers installed on the Planck spacecraft. The pointing for each 
bolometer is computed by combining the attitude with the loca- 
tion of the bolometer in the focal plane reconstructed from Mars 
observations. 

From the cleaned TOI and the pointing, channel maps have 
been made using bolometers at a given frequency. The path from 
TOI to maps in the HFI DPC is schematically divided into three 
steps, ring-making, destriping and map-making. The first step 
averages circles within a pointing period to make rings with 
higher signal-to-noise ratio taking advantage of the redundancy 
of observations provided by the Planck scanning strategy. The 
low amplitude 1 // component is accounted for in a second step 
using a destriping technique. Finally, cleaned maps are produced 
using a simple co-addition of the rings. 



The noise in the channel maps is essentially white with 
a mean standard deviation of 1.4 x 10~ 3 , 4.1 x 10~ 3 , 1.4 x 
10 3 MJy/sr at 353, 545 and 857 GHz respectively (Planck HFI 
Core Team 201 lb). The photometric calibration is performed ei- 
ther at the ring level using the CMB dipole, for the lower fre- 
quency channels, or at the map level using FIRAS data, for the 
higher frequency channels at 545 and 857 GHz. The absolute 
gain calibration of HFI Planck maps is known to better than 2% 
at 353 GHz and 7% at 545 and 857 GHz (see Table 2 in Planck 
HFI Core Team 2011b). 

The detection algorithm requires the use of ancillary data to 
trace the warm component of the gas. Thus we combine Planck 
data with the IRIS all-sky data (Miville-Deschenes & Lagache 
2005). The choice of the IRIS 100 pm as the warm template is 
motivated by the following: (i) 100 pm is very close to the peak 
frequency of a black body at 20 K, and traces the warm compo- 
nent of the Galaxy; (ii) the fraction of small grains at this wave- 
length remains very small and does not significantly the estimate 
of the emission from large grains that is extrapolated to longer 
wavelengths; (iii) the IRAS survey covers almost the entire sky 
(only 2 bands of ~2% of the whole sky are missing); (iv) the 
resolution of the IRIS maps is similar to the resolution of Planck 
in the high frequency bands, i.e. around 4.5'. Using the map at 
100 pm as the warm template is, of course, not perfect, because 
a non-negligible fraction of the cold emission is still present at 
this frequency. This lowers the intensity in the Planck bands af- 
ter removal of the extrapolated background. We will describe in 
detail, especially in Sect. 2.3, how we deal with this issue for the 
photometry of the detected cores. 

All Planck and IRIS maps have been smoothed at the same 
resolution 4.5' before source extraction and photometry process- 
ing. 

2.2. Source Extraction Method 

We have applied the detection method described in Montier et 
al. 2010, known as CoCoCoDeT (standing for Cold Core Colour 
Detection Tool), on the combined IRIS plus Planck data set de- 
scribed in Sect. 2.1. This algorithm uses the colour properties of 
the objects to be detected to separate them from the background. 
In the case of cold cores, the method selects compact sources 
colder than the surrounding envelope and the diffuse Galactic 
background, that is at about 17 K (Boulanger et al. 1996) but 
can largely vary from one place to the other across the Galactic 
plane or at higher latitudes. This Warm Background Subtraction 
method is applied on each one of the three Planck maps, and 
consists of 6 steps: 

1 . for each pixel, the background colour is estimated as the me- 
dian value of the Planck map divided by the 100pm map 
within a disc of radius 15' around the central pixel; 

2. the warm component in a pixel at the Planck frequency is ob- 
tained by multiplying the estimate of the background colour 
with the value of the pixel in the 100pm map; 

3. the cold residual map is computed by subtracting the warm 
component from the Planck map; 

4. the local standard deviation around each pixel in the cold 
residual map is estimated in a radius of 30' using the so- 
called Median Absolute Deviation that ensures robustness 
against a high confusion level of the background and pres- 
ence of other point sources within the same area; 

5. a thresholding detection method is applied in the cold resid- 
ual map to detect sources at a signal-to-noise ratio SNR>4; 
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6. final detections are defined as local maxima of the SNR con- 
strained so that there is a minimum distance of 5' between 
them. 

This process is performed at each Planck band yielding in- 
dividual catalogues at 857 GHz, 545 GHz and 353 GHz. The 
last step of the source extraction consists in merging these three 
independent catalogues requiring a detection in all three bands 
at SNR>4. This step rejects spurious detections that are due to 
map artifacts associated with a single frequency (e.g. stripes or 
under-sampled features). It increases the robustness of the final 
catalogue, which contains 10783 objects. 

We stress that no any other a-priori constraints are imposed 
on the size of the expected sources, other than the limited area 
on which the background colour is estimated. Thus the maxi- 
mum scale of the C3PO objects is about 12'. Note also that this 
Warm Background Subtraction method uses local estimates of 
the colour, identifying a relative rather than an absolute colour 
excess. Thus cold condensations embedded in cold regions can 
be missed, while in hot regions condensations may be detected 
that are not actually cold. A more detailed analysis in tempera- 
ture is required to assess the nature of the objects. 

2.3. Photometry 

We have developed a dedicated algorithm to derive the photom- 
etry of the clump itself. The fluxes are estimated from the cold 
residual maps, instead of working on the initial maps where the 
clumps are embedded in their warm surrounding envelope. As 
already stressed above, the main issue is to perform the photom- 
etry on the IRIS 100 pm maps that also include a fraction of the 
cold emission. The flux of the source at 100 pm has to be well 
determined for two reasons: (1) an accurate estimate of the flux 
at this frequency is required because it is constrains significantly 
the rest of the analysis (in terms of spectral density distribution 
(SED) and temperature); (2) an incorrect estimate of the flux at 
100 pm will propagate through the Planck bands after removal 
of the extrapolated warm component. The main steps of the pho- 
tometry processing are described in the following subsections. 
An illustration of this process is provided in Fig. B.5 of the asso- 
ciated Planck Early Paper on Cold Clumps describing in detail a 
sample of 10 sources (Planck Collaboration 201 lr). 

2.3.1 . Stepl : Elliptical Gaussian fit 

An elliptical Gaussian fit is performed on the 1° x 1° colour map 
857 GHz divided by 100 pm centered on each C3PO object. This 
results in estimates of three parameters: major axis extension 
0"Maj, minor axis extension cr Min and position angle if/. The re- 
lation between the extension cr and the FWHM 9 of a Gaussian 
is given by : 

a = 6/ V81n(2) (1) 

If the elliptical Gaussian fit is indeterminate, a symmetrical 
Gaussian is assumed with a FWHM fixed to 6 = 4.5', and the 
flag Aper Forced is set to on. In these cases, the source fluxes are 
severely underestimated at all frequencies. This flagged popula- 
tion contains 978 sources which are rejected from the physical 
analysis of Sect. 4, but not from the entire catalogue, which is 
used to assess the association with ancillary data (cf Sect. 3) and 
to study morphology at large scale (cf Sect. 5). 



2.3.2. Step2: 100yum photometry 

The photometry on the 100 pm map is obtained by surface fit- 
ting, performed on local maps of 1° x 1° centered on each can- 
didate. All components of the map are fitted as a whole: a poly- 
nomial surface of an order between three and six for the back- 
ground; a set of elliptical Gaussians when other point sources are 
detected inside the local map; and a central elliptical Gaussian 
corresponding to the cold core candidate for which the elliptical 
shape is set by the parameters obtained during step 1 . When the 
fit of the background is poor, i.e. a clear degeneracy is observed 
between the polynomial fit and the central Gaussian, we switch 
to a simple aperture photometry on the local map. Note that the 
aperture photometry is performed taking into account the ellipti- 
cal shape of the cold core provided by step 1. In such cases (140 
sources), the flag Bad Sfit 100pm is set to on. Occasionally no 
counterpart at all is observed at 100 pm, when the cold core can- 
didate is too faint or very cold, or the confusion of the Galactic 
background is too high. In such case, we are not able to derive 
any reliable estimate of the 100 pm flux of the core, so only an 
upper-limit can be provided. This upper limit is defined as three 
times the standard deviation of the cold residual map within a 
25' radius circle, and the flag Upper 100pm is set to on. There 
are 2356 objects for which only an upper limit of the temper- 
ature is derived. This population represents a very interesting 
sub-sample of the whole catalogue, probably the coldest objects, 
but we do not have confidence in the physical properties derived 
from the Planck data and so it is excluded from the physical 
analysis. 

2.3.3. Step 3: lOO^m correction 

Once an estimate of the flux at 100 pm has been provided by 
steps 1 and 2, the warm template at 100 pm is corrected by re- 
moving an elliptical Gaussian corresponding to the flux of the 
central clump. This new warm template is then extrapolated and 
subtracted from the Planck maps to build the cold residual maps. 
When only an upper limit has been obtained at 100 pm, the warm 
template is not changed. 

2.3.4. Step 4: Planck bands photometry 

Aperture photometry is performed on local cold residual maps 
centered on each candidate in the Planck bands, at 857 GHz, 
545 GHz and 353 GHz. This aperture photometry takes into ac- 
count the real extension of each object by integrating the sig- 
nal inside the elliptical Gaussian constrained by the parameters 
obtained at step 1. The background is estimated by taking the 
median value on an annulus around the source. Nevertheless, in 
229 cases, no positive estimate of the flux has been obtained, be- 
cause of the presence of cold point sources that are too close or 
because the background is highly confused. These sources (for 
which the flag PS Neg is set to on) are simply removed from the 
physical analysis described in this paper. 

2.4. Monte-Carlo Quality Assessment 

To assess the quality of our photometry algorithm, we have 
performed a Monte-Carlo analysis. A total of 10000 simulated 
sources are randomly distributed over the whole sky in the IRIS 
and Planck maps. The sources are assumed to follow the emis- 
sion of a modified black body with a temperature randomly, T, 
distributed between 6 K and 20 K, and an associated spectral in- 
dex given by B = 11.5 x T-° bb within a 20% error bar, based 
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lcr(%) 


Bias (vc) 


lcr(%) 


Flux at 100 /im 


1.4 


31.7 


1.0 


4.7 


-58.1 


14.1 


i n i 
1 1 /. 1 


1 OA A 


Flux at 857 GHz 


-5.0 


6.2 


3.9 


3.2 


-56.3 


13.8 


-11.0 


6.0 


riuA at jh-j v j i 1 / 


-3.6 


6.4 


3.7 


3.7 


-55.8 


14.4 


q n 

y.\j 


ft 


Flux at 353 GHz 


-5.0 


7.3 


2.4 


4.7 


-58.9 


14.8 


-10.0 


6.7 


FWHM 


-0.6 


16.2 


30.9 


27.7 


-25.2 


16.3 


-6.7 


15.3 


Ellipticity 


0.0 


8.2 


0.0 


9.5 






0.0 


9.0 


T 


-4.2 


5.2 


-4.1 


1.6 


-6.5 


3.8 


0.4 


16.0 


P 


9.8 


7.3 


10.5 


2.4 


11.2 


6.7 


2.7 


18.7 



Table 1. Statistics of the Monte-Carlo analysis performed to estimate the robustness of the photometry algorithm. The bias (ex- 
pressed in %) is defined as the relative error between the median of the output distribution of the photometry algorithm and the 
injected input. The l<x (expressed in %) represents the discrepancy around the most probable value of the output distribution. Those 
quantities are given in the various cases corresponding to the output flags provided by the algorithm. Statistics of the temperature 
and spectral index is also given here to show the impact of the observed error on fluxes. 





Normal 


Bad Sfit 100 pm 


Aper Forced 


Upper 100pm 


Quantity 


Bias (%) 


lcr{%) 


Bias(%) 


lo-(%) 


Bias(%) 


lcr(%) 


Bias (%) 


lcr(%) 


Flux at 100 yum 


11.5 


44.3 


0.8 


8.4 


-51.6 


21.1 


204.5 


278.2 


Flux at 857 GHz 


-4.0 


8.1 


2.1 


4.7 


-58.3 


20.1 


-10.4 


7.1 


Flux at 545 GHz 


-2.5 


8.0 


2.4 


4.9 


-57.4 


21.3 


-7.8 


7.0 


Flux at 353 GHz 


-3.4 


8.7 


1.9 


5.5 


-59.3 


21.3 


-8.7 


7.4 


FWHM 


0.0 


18.1 


31.0 


31.1 


-24.4 


16.9 


-5.2 


17.6 


Ellipticity 


0.0 


9.3 


-0.5 


9.2 






0.1 


10.4 


T 


-2.1 


6.3 


-3.2 


1.8 


-4.4 


6.2 


6.8 


20.6 


P 


7.1 


8.2 


9.3 


2.6 


5.6 


12.3 


-4.9 


20.4 



Table 2. Same as Table 1 in the Galactic plane (|b| < 25°). 



on the work done on Archeops data by Desert et al. (2008). The 
FHWM of the simulated sources spans from 4.5' to 7' with an 
ellipticity ranging from to 0.87. The flux at 857 GHz is taken 
from 10 to 500 Jy following a logarithmic random distribution. 
The derived fluxes in all IRAS and Planck bands take into ac- 
count the colour correction. We apply our complete process of 
photometry on this set of simulated data, and retrieve an esti- 
mate of all quantities (fluxes, FWHM, ellipticity) in the various 
cases described by the flags listed before (cf Fig. A. 1). Statistical 
bias and Icr errors are derived for all quantities and cases, and 
are listed in Table 1 and 2 for all-sky and \b\ < 25° respectively. 
We the temperature and spectral index estimates recovered at the 
end of the processing are also listed to illustrate the impact of the 
errors on the fluxes. 

This Monte-Carlo analysis confirms, firstly, why sources 
with Aper forced set to on should be rejected from the physical 
study, since for these sources fluxes are systematically under- 
estimated by about 60%. Sources with Upper 100 pm set to on, 
for which only an upper limit at 100 pm has been provided by 
the algorithm, the flux at 100 pm is over-estimated by a factor 
of two, with an associated discrepancy that can reach a factor of 
three times the input value in regions close to the Galactic plane. 
Moreover the fluxes in the Planck bands are significantly biased 
to lower values, with a bias greater than the Icr discrepancy. The 
resulting temperature estimate is, as expected, greater than the 
injected value and the uncertainties in the temperature and spec- 
tral index are around 20%. This illustrates the limitations on any 
physical conclusions that could be drawn from this population 
of sources. When a bad fit of the 100 pm background has been 
obtained, Bad Sfit 100 pm flag set to on, the main error comes 
from the highly biased estimate of the FWHM (~31%), leading 
to an over-estimate of the fluxes in all bands. This happens when 



a strong source is embedded in a faint background (e.g. at high 
latitude), introducing a degeneracy between the fit of the central 
elliptical Gaussian and the polynomial fit of the background sur- 
face at 100 pm. Although bias and l-cr values are smaller than 
in the normal case due to the strong signal of these sources, we 
reject this population from the physical analysis, because they 
could introduce wrong estimates of the physical properties based 
on a highly biased extension. 

If we focus now on the normal case, when the photometry 
algorithm has performed well, we first observe a slight bias of 
all fluxes estimates. The bias at 100 pm becomes larger when 
looking into the Galactic plane (11.5% for |b| < 25 compared 
to 1.4% over the whole sky). The fluxes Planck bands, however, 
are less under-estimated when looking inside the Galactic plane, 
with biases spanning from 2.5% to 5%. The associated Icr er- 
rors are about 6 to 7% on all-sky and 8-9% in the Galactic plane. 
The impact of such a biased estimate of the fluxes will be dis- 
cussed together with the study on the calibration uncertainty in 
Sect. 4. 1 . On the other hand, the FWHM estimate are typically 
biased by less than 1% and have an accuracy of ~18%, when the 
ellipticity presents no bias and an accuracy of ~9%. Finally the 
temperature and spectral index are derived using the method de- 
scribed in Sect. 4. 1 . Whereas the temperature is slightly under- 
estimated (~2% in the Galactic plane), the associated spectral 
index is over-estimated by ~7%. The statistical l-cr uncertain- 
ties are about 6% and 8% for T and B respectively. These results 
will be taken into account in detail when discussing the physical 
properties of these cold sources in Sect. 4. 1 . 

The Monte-Carlo simulations described here demonstrate 
the robustness of our photometry algorithm, and justify the re- 
jection of entire categories of objects using the photometry flags, 
such as the Aper Forced, PS Neg and Bad Sfit 100pm. The re- 
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Fig. 1. Colour-Colour diagram of the catalogue. The over-plotted 
symbols stand for the positive cross-matches with non ISM ob- 
jects. The red contours give the domain of the diagram filled by 
Archeops cold cores assumed to follow a grey-body law, with a 
temperature ranging from 6 K < T < 25 K, and a spectral index 
P given by Desert et al. (2008). 

maining sample consists of 9465 objects, divided into two cat- 
egories: 1840 objects have only an upper limit estimate of the 
flux at 100 /urn and 7625 have well defined photometry in IRAS 
and Planck bands. We will focus on this last category of 7625 
sources for the rest of the analysis on the physical properties. 
Based on this Monte-Carlo analysis, we will adopt the following 
estimate of the lcr uncertainty on fluxes: 40% on IRAS 100 fim, 
and 8% on Planck bands. This error is much larger than the in- 
trinsic pixel noise and so instrumental errors are neglected. 



1200 




SNR 

Fig. 2. Signal-to-noise ratio (SNR) of new sources (dash line) 
overlaid on the SNR of all sources (solid line). 



Simbad type 


C3PO 


<MC > 




[%] 


[%] 


ISM 


49.0 


21.7 


Star 


2.3 


4.9 


Gal 


2.1 


7.4 


Radio 


5.3 


7.7 


QSO 


0.1 


0.3 


Others 


0.3 


0.2 


New detections 


40.9 


57.8 



Table 3. Cross match with Simbad database for C3PO and simu- 
lated catalogues, for each category of Simbad type. The < MC > 
column gives an estimate of the probability of chance alignment 
for each Simbad type. 



2.5. Cross-Correlation with existing catalogues 

As one step of the validation of our detections, we have per- 
formed an astrometric search on the Simbad database 2 for all 
known sources within a 5' radius of our sources. There are a 
large number of objects in the Simbad database which raises 
the question of chance alignments. This is especially true for 
extragalactic objects which have a reasonably isotropic sky dis- 
tribution. To judge the number of chance alignments that can 
be expected by performing this kind of search, we have also 
conducted a Simbad cross check on the positions of a set of 
100 Monte-Carlo simulated catalogues presented in Sect. 5.1.1. 
These Monte-Carlo realizations reproduce the object density of 
the Planck catalogue per bin of longitude and latitude. The re- 
sults presented in Table 3 show that the number of coincidences 
in the ISM category is greater in the C3PO catalogue than the 
probability of chance alignment estimated from the Monte-Carlo 
simulations. On the contrary, the fraction of contaminants (i.e. 
Galaxies, QSO, Radio Sources, stars) is always lower in C3PO 
than in the Monte-Carlo realizations. Thus extragalactic objects 
and Galactic non-dusty objects are mostly rejected by the de- 
tection algorithm, whereas actual ISM structures are preferen- 
tially detected. A more detailed comparison between C3PO and 
IRDCs catalogues is presented in Sect. 7.1. 

Nevertheless the association with probable contaminants in 
C3PO is quite high (~10%) and not all are necessarily the result 
of chance alignments. To disentangle between chance alignment 
and real matches, we use colour-colour information as shown 
in Fig. 1. Mostly objects are distributed in the bottom-left cor- 
ner of the diagram, typical of dust-dominated emitters. The red 

2 http://simbad.u-strasbg.fr/simbad/ 



contours of this figure show the domain filled by dusty objects 
assuming a grey-body emission law, with 6 K < T < 25 K, and 
a spectral index /3 given by Desert et al. (2008). The match be- 
tween Planck detections and this colour-colour domain is strong. 
Only a few objects (17) show the colour-colour properties of 
radio emitters, located in the top-right corner, indicating real 
matches with extragalactic objects. For the rest of the sample, 
the probability of chance alignment is high. Concerning the as- 
sociation with stars, except for a few X-ray emitters, mostly all 
Simbad matches seem associated with dusty emission, and thus 
represent chance alignment. 

We finally reject only the obvious extragalactic matches, lo- 
cated in the top-right corner of the colour-colour diagram, lead- 
ing to 7608 objects. 

Out of the 7608 sources in the photometric reliable cata- 
logue, 40 % have no counterpart in the Simbad database. In ad- 
dition, these new detections have a similar SNR distribution as 
the entire catalogue as shown in Fig. 2, and can be considered as 
reliable as the entire catalogue. 

3. Spatial Distribution 

3.1. Association with Galactic structures 

The all-sky distribution of the 10783 C3PO sources is presented 
in the upper panel of Fig. 3. Mostly concentrated in the Galactic 
plane, the distribution clearly follows Galactic structures be- 
tween latitudes of -20° and +20°. A few detections are observed 
at high Galactic latitude (|b| > 30°) and after cross-correlation 
with external catalogues have been confirmed not to be extra- 
galactic objects (see Sect. 2.5). 
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Fig. 3. Upper panel: All-sky map of the number of C3PO Planck cold clumps per sky area, smoothed at 3°. Middle panel: CO 
contours are over-plotted on the C3PO density map which is set to where CO map is not defined. Lower panel: Av contours are 
over-plotted on the C3PO density map which is set to where Av map is lower than 0. 1 Av. 
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In the middle panel of Fig. 3, contours of the integrated in- 
tensity map of the CO J 1-0 line are overlaid on the Planck cold 
clumps density all-sky map. This CO map is a combination of 
CO data from Dame et al. (2001 ) and NANTEN data (Fukui et al. 
1999; Matsunaga et al. 2001; Mizuno & Fukui 2004), as defined 
in Planck Collaboration (201 lo). The correlation between CO 
and C3PO Cold Clumps is quite impressive and demonstrates 
once again the robustness of the detection process and the con- 
sistency of the physical nature of these Planck cold objects. A 
detailed analysis shows that more than 95% of the clumps are 
associated with CO structures. 

The lower panel shows the same kind of spatial correlation 
with the all-sky Av map (Dobashi 201 1 in preparation). The Av 
map traces more diffuse regions of the Galaxy and extents to 
higher latitude, where cold clumps are also present. About 75% 
of the C3PO objects are associated with an Av signature greater 
than 1. 



Name 


Lon 


Lat 


Area 


Distance 


Nb 




[deg] 


[deg] 


[deg 2 ] 


[pc] 




Aquila Serpens 


3 


28 


30 


260 


59 


Polaris Flare 


24 


123 


134 


150 


55 


Camelopardalis 


20 


148 


159 


240 


11 


Ursa Major 


35 


148 


44 


240 


13 


Taurus 


-15 


170 


883 


140 


393 


Taurus Perseus 


-15 


170 


883 


350 


227 


A Ori 


-13 


196 


113 


400 


66 


Orion 


-9 


212 


443 


450 


353 


Chamaeleon 


-16 


300 


27 


150 


114 


Ophiuchus 


17 


355 


422 


150 


311 


Hercules 


9 


45 


35 


300 


16 



Table 4. Molecular complexes used to associate C3PO cold 
clumps to Galactic well-known structures, for which an estimate 
of the distance is available. 



3.2. Distance Estimation 

Distance estimates are essential to properly analyse the popula- 
tion of detected cold clumps. We have used four different meth- 
ods: association with IRDCs, association with known molecular 
complexes, three dimensional extinction method using 2MASS 
data, and extinction method using SDSS data. 

3.2.1. Distances to IRDCs 

Simon et al. (2006b) and Jackson et al. (2008) provide kinematic 
distance estimates for a total of 497 IRDCs extracted from the 
MSX catalogue (Simon et al. 2006a) that consists of 10931 ob- 
jects. Kinematic distances are obtained via the observed radial 
velocity of gas tracers in the plane of the Galaxy. By assum- 
ing that the Galactic gas follows circular orbits and a Galactic 
rotation curve, an observed radial velocity at a given longitude 
corresponds to a unique Galactocentric radius. Of course, this 
means that in the inner Galaxy, two heliocentric distances are 
possible. This technique is only applicable in the plane and re- 
quires the availability of appropriate molecular data. We find 
127 Planck cold clumps, over the complete catalogue, associ- 
ated with IRDCs that already have a kinematic distance estimate. 
This number decreases to 32 associations over the 7608 objects 
of the photometric reliable C3PO catalogue. 

A more recent work by Marshall et al. (2009) uses an extinc- 
tion method, detailed in Sect. 3.2.3, on the same MSX catalogue 
of IRDCs to derive the distance of 1259 objects. This yields 188 
associations with C3PO clumps over the entire catalogue, and 
47 over the photometric reliable C3PO catalogue. 

3.2.2. Distances to known molecular complexes 

The all-sky distribution of cold clumps follows known molec- 
ular complexes. Many of these have distances estimates in the 
litterature. To assign the distance of a complex to a particular 
cold clump we use the CO map of Dame et al. (2001) to trace 
the structure of the molecular cloud above a given threshold, and 
test for the presence of cold clumps inside this region. The as- 
sociation has been performed on 14 molecular complexes (see 
Table 4), leading to 1152 distance estimates over the entire cat- 
alogue and 947 on the photometrically reliable catalogue, cata- 
logue. 



3.2.3. Distances from extinction signature 

Genetic forward modelling (using the PIKAIA code 
Charbonneau 1995) is used along with the Two Micron 
All Sky Survey (Skrutskie et al. 2006) and the Besancon 
Galactic model (Robin et al. 2003) to deduce the three di- 
mensional distribution of interstellar extinction towards the 
cold clump detections. The derived dust distribution can then 
be used to determine the distance and mass of the sources, 
independently of kinematic models of the Milky Way. Along 
a line of sight that crosses a cold clump, the extinction is seen 
to rise sharply at the distance of the cloud. The method is fully 
explained in Marshall et al. (2006) and Marshall et al. (2009). 

The distance, as determined by this technique, provides line 
of sight information on the dust distribution. However, it does 
not have sufficient angular resolution to perform morphologi- 
cal matches on the cold clumps. To ensure that the extinction 
rise detected along the line of sight is indeed related to the inner 
structure we perform a consistency check on the column density 
derived from the extinction and from the source flux, corrected 
for its temperature. Only detections where the two column den- 
sities are in agreement within a factor of two are retained. This 
leads to distance estimates for 978 objects of the entire and pho- 
tometric reliable catalogue. 

3.2.4. Distances from SDSS 

Distances to cold clumps within 1 kpc are obtained by analysis 
of distance-reddening relations for late spectral type stars within 
the line of sight to each source (Mc Gehee 201 1 in preparation). 
Specifically, we use Sloan Digital Sky Survey photometry of Ml 
to M5 dwarfs colour-selected by the reddening-invariant index 

Qgri = (g-r)- E _ . (r - i). (2) 

The updated ugriz reddening coefficients of Schlafly et al. (2010) 
are used. The median stellar locus of Covey et al. (2007) forms 
the basis of a calibration between Q gr j and the intrinsic g — i 
colour. After dereddening, the distance to each star is determined 
including corrections for Galactic metallicity variation following 
Bochanski et al. (2010). 

The distance-reddening profile is constructed by computing 
the median reddening for stars within a circular patch centered 
on the core location for 25 pc wide distance bins spanning 
to 2000 pc. We fit the observed reddening profile to the model 
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Fig. 4. Distribution of C3PO cold clumps as seen from the North 
Galactic Pole. Colours stand for methods used to estimate dis- 
tance: Molecular Complex association (green), SDSS extinc- 
tion (light blue), 2MASS extinction (dark blue), IRDCs extinc- 
tion (orange) and IRDCs kinematic (red). The red dashed circle 
shows the 1 kpc radius around the sun. Black dashed lines rep- 
resent the spiral arms and local bar. The black circles give the 
limits of the molecular ring. 

defined by convolution of the near-field plus single cloud profile 
with a Gaussian (in distance modulus), this function is: 



E(B - V). 



obs 



a + c 



f 



l 



(3) 



where x is the independent variable (distance modulus), xq is the 
location of the single cloud, a is the near-field reddening, c is the 
reddening associated with the cloud, and cr is the width of the 
Gaussian. The fitted cr values are typically 0.4 to 0.5 magnitudes 
in m — M, as expected from the standard deviation of the (r - 
z, M r ) used to assign absolute magnitudes. 

Analysis of calibration fields containing well-studied molec- 
ular cloud complexes, e.g. the Orion B Cloud, reveal that the 
recovered distance moduli are underestimated by 0.2 to 0.3 mag- 
nitudes, consistent with the bias expected from the M dwarf mul- 
tiplicity fraction. 

This processing leads to 1452 distance estimates over the en- 
tire catalogue and 1004 over the photometric reliable one. 

3.2.5. Combined results 

The number of sources for which distances could be recovered 
depends on the method used (cf Table 5). There is some over- 
lap but each method has its distinct advantages according to the 
distance range being considered. The 2MASS extinction method 
is not very sensitive nearby (D<1 kpc), as there are not enough 
stars to determine accurately the line of sight information. In 
contrast, the extinction method using SDSS is especially de- 
signed for nearby objects. For objects with 1 kpc, we have used 
SDSS distances when availble or molecular complex distances. 



Method 


hntire CjrU 


Reduced LorU 




/ 1 (Y7 OQ\ 

(lU/oo) 




1 IX U\^, N ^JVlIlClllaLlL,^ 


1 97 
1 _ / 


jz 


IRDCs (Extinction) 


188 


47 


2MASS Extinction 


978 


978 


SDSS Extinction 


1452 


1004 


Molecular Complexes 


1152 


947 


Total 


3411 


2619 



Table 5. Number of distance estimates available of the C3PO 
sources for each method. Notice that the total numbers are not 
equal to the sum of all methods, due to overlap between them. 
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Fig. 5. Distance distribution of the MSX IRDCs (Simon et al. 
2006b) (solid line) and of the subset associated to the cold 
clumps of the entire C3PO catalogue (dot-dash-dash line) and 
the photometric reliable subset of C3PO (dotted line). 



The number of objects for which we have a distance estimate 
is 2619 out of a total of 7608 objects in our photometric reli- 
able subset, i.e. ~34%. The distances of the cold clumps span 
from 0.1 to 7 kpc, but they mainly concentrated in the nearby 
Solar neighbourhood as shown on Fig. 4. This type of distri- 
bution has been already demonstrated using simulations, see 
Fig. 10 of Montier et al. (2010). The lack of detections at large 
distances is mainly caused by the effects of confusion within 
the Galactic plane, from which suffers the detection method. 
Nevertheless, when comparing the distance distribution of the 
C3PO cold clumps associated to MSX IRDCs with the total sam- 
ple of Simon et al. (2006b) in Fig. 5, we notice that the fraction 
of C3PO - IRDCs matches does not depend on distance and ex- 
tends to 8 kpc. 

Because the subset of C3PO cold clumps with a distance es- 
timate has been obtained using different methods, exploring var- 
ious regions and distances over the sky, this sample appears het- 
erogeneous. The completeness of the catalogue with distances is 
quite difficult to assess. Thus we define two subsets for further 
analysis, especially when looking at number counts, for which 
we know that the sample is more homogeneous: the first subset 
(1790 objects) deals with the local objects (D < 1 kpc) and uses 
only estimates from molecular complexes association and SDSS 
extinction; the second subset (674 objects) focuses on distant ob- 
jects (D > 1 kpc) and uses only 2MASS extinction estimates and 
IRDCs associations. 
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Fig. 6. Distribution of the temperature of the cold clumps (blue), 
of the warm envelope (red) and of the total (green) estimated 
inside the elliptical Gaussian of the clump itself. The averaged 
temperature of the local background is plotted in red dot-dash 
line. 



4. Physical Properties 

4.1. Temperature 

The temperature of the sources is estimated from SEDs using 4 
bands: the IRAS 100 /urn and the three highest frequency Planck 
bands 857 GHz, 545 GHz and 353 GHz. The assumed emission 
model is a modified black-body law, defined as: 



S v =AB v (T)v p , 



(4) 



where S v is the flux integrated over the solid angle 
Clc = TrcMajCMin, A is the amplitude, T is the temperature, 
B is the spectral index and B v is the Planck function. 

For each source, a set of four temperatures is measured: (1) 
the temperature of the clump Tq is defined as the temperature 
based on the SEDs of the cold residual as described in Sect. 2.3; 
(2) the temperature of the warm envelope T env is obtained from 
aperture photometry over the same region but performed on the 
warm component; (3) the total temperature T tot is defined as the 
temperature of the source in the initial map, i.e. without remov- 
ing any warm component; (4) the temperature of the local back- 
ground Tbkg is defined as the temperature of the average surface 
brightness around the source. 

We have first fixed the spectral index to B — 2 (Boulanger 
et al. 1996). Ax 2 fit is performed on the SEDs to derive all es- 
timates of temperatures and associated 1-cr errors. The distribu- 
tion of these temperatures is shown on the upper panel of Fig. 6. 
The temperature of the cores Tc (blue line) peaks at 13.4K and 
spans from 9 K to 16 K. The temperature of the total T, , (green 
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Fig. 7. Reduced x 2 obtained in the case B — 2 as a function of 
the temperature obtained with B free. When T becomes lower, 
the x 2 becomes larger. 



line), of the warm envelope T env (red line), and of the local back- 
ground T^g (red dot dash line) distributions peak respectively at 
13.9K, 15.1 K and 16.1 K. The uncertainty on the temperature 
estimates is about 7%. These results are in good agreement with 
the expected values of cold cores (e.g. Bergin & Tafalla 2007) 
and consistent with the results of our Monte-Carlo simulations 
demonstrating that our source extraction method accurately re- 
covers the cold source parameters in the presence of a warmer 
background. 

In a second analysis, we performed a three parameter (A, T 
and B) x 2 fit leading to the temperature distributions shown in 
the lower panel of Fig. 6. The^ 2 fit is performed on a grid tak- 
ing into account the colour correction as defined in Planck HFI 
Core Team (2011b) and gives the exact minimum of the x 2 in 
the (A,T, B) space and providing the associated 1-cr uncertainty. 
Evaluating the^f 2 obtained with B = 2 as a function of the best fit 
temperature obtained from the full three parameter fits, we see 
that a model B — 2 is reasonable for temperatures in the range 
10K<T<18K (for which the^ 2 < 1), but does not provide 
a good fit at lower temperature T < 10 K (see Fig. 7). In fact, 
the lower the temperature, the worse the fit. Using B as a free 
parameter, the temperature distributions peak at 13 K, 13.9 K, 
15.5 K and 17 K for Tc, T tot , T env and Tbkg respectively, with an 
error of about 7%. The associated spectral index B varies from 
1.5 to 3, with an uncertainty of 21% and a mean value of 2.1 for 
cold clumps and 1.8 for the total emission, consistent with other 
studies based on Planck data (Planck Collaboration 2011o,t,u). 
The temperature of the cold clumps span the range 7 K to 17 K. 

The bias and the uncertainty of the temperature and spectral 
index have to be adjusted, taking into account the Monte-Carlo 
analysis of the photometry algorithm (see Sect. 2.4), and the im- 
pact of the calibration uncertainty detailed in Sect. B. We recall 
that a bias of — 2% on T and ~7% on B is induced by the pho- 
tometry itself. On the other hand, the calibration uncertainty of 
fluxes does not introduced any bias on T or B, but generates an 
error of ~8% onB and from 3% to 5% on T, that should be added 
quadratically to the uncertainty due to statistical errors. All these 
considerations lead to a final range of temperature spanning the 
range 7 K to 17 K with an uncertainty of about 9%, and a spec- 
tral index B varying from 1.4 to 2.8 with an uncertainty of about 
23%. 
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Fig. 8. Upper panel: distribution of the FWHM of the Planck de- 
tections compared to the distribution of the local PSF at 857 GHz 
(dashed line). Middle panel: distribution of the ellipticity of the 
cold clumps (solid line) and of the local PSF (dashed line). 
Lower panel: distribution of the position angle of the elliptical 
Gaussian of the clumps (solid line), of the local PSF (dashed 
line), and difference between both (dotted line). 

4.2. Extension and ellipticity 

The extension and ellipticity of the sources derived during step 
1 of the photometry algorithm described in Sect. 2.3 have been 
compared with the local Point Spread Function (hereafter PSF) 
provided by the FEBeCoP tool (Mitra et al. 2010) at 857 GHz. 
This PSF takes into account the scanning strategy and the pix- 
elization of the maps at each location of the sky. Thus for each 
source, an elliptical Gaussian fit is applied on the PSF smoothed 
at 4.5' to get the local FWHM 6>psf, ellipticity epsF and position 
angle i^psf of the effective beam. The FWHM 8 is defined as the 
geometric mean of the major and minor axis widths: 



8 - -^(#Maj ' #Min), 

and the ellipticity is given by: 



e= Jl 



?Min 



(5) 



(6) 



A few examples of FEBeCoP beams for Planck HFI detec- 
tors are given in Fig. B.5 of Planck Collaboration (201 lr). 
Fig. 8 compares the statistical distributions of the FWHM (up- 
per panel), ellipticity (middle) and position angle (lower panel) 
between C3PO sources (solid line) and the local PSF at 857 GHz 
(dashed line). 



Cold clumps are clearly extended, with an average value of 
8c of 7.7' compared to the 4.3' of the average PSF over the sky. 
Assuming that these compact sources are resolved by the Planck 
beam, we can deconvolve them to derive the inferred intrinsic 
source size (dot-dash line in Fig. 8): 



4' 



(fi 

°PSF' 



(7) 



where 8 C is the extension of the source and 8 PS F is the PSF ex- 
tension. We find that 8j/8psF ~ 1-4, and could conclude that we 
have resolved the sources. Nevertheless, as pointed out by Enoch 
et al. (2007) and Netterfield et al. (2009), the fact that the cold 
sources are mostly extended compared to the PSF is an indica- 
tor of the hierarchical structure of these objects. Netterfield et al. 
(2009) show that the BLAST sources present the same behav- 
ior with a ratio between the inferred source size and the BLAST 
beam equal to 1.1. Enoch et al. (2007) obtained a value of 1 .5 
for cold cores in Serpens, Perseus and Ophiucus observed with 
Bolocam. Indeed these compact sources are associated to larger 
envelopes presenting radial density profiles in power law with an 
exponent equal to -2 to -1 (Young et al. 2003). 

Cold clumps are also mostly elongated, with a distribution 
of axial ratios extending to values as large as 5 and peaking at 
around 1.5, compared to the mean value of 1.3 for the local PSF. 
Note that the C3PO cold clumps are not preferentially aligned 
with the major axis of the PSF, the position angles of the ellip- 
tical clumps and of the PSF are uncorrelated. As also stressed 
by Planck Collaboration (201 lr), cold cores are often associ- 
ated with filaments and parts of larger elongated cold structures 
where star formation occurs. This was noted a long time ago by 
Barnard (1907) for Taurus, and it has more recently been inves- 
tigated by Herschel observations in Polaris (Men'shchikov et al. 
2010) and Aquila (Konyves et al. 2010). This characteristic of 
the cold core population can now addressed more generally us- 
ing the Planck all-sky data and is discussed in detail in Sect. 5. 

When distances are available (see Sect. 3.2), we can derive 
the physical size of the sources, defined as the FWHM in pc. 
Fig. 9 presents the statistical distribution of the size obtained for 
2619 sources. A distinction is made between local (D < lkpc, 
dashed line ) and distant (D > lkpc, dot-dash line) sources as 
defined in Sect. 3.2.5. The relation of Elmegreen & Falgarone 
( 1996) and Heyer et al. (2001), a size spectral index of -2.3 typi- 
cal of dust clouds, is over-plotted on the distributions of the two 
subsets. 

4.3. Column Densities 

The column density values averaged over the clump solid angle 
can be derived from the integrated flux using : 



Qc/Uf«HKv x B Vo (Ty 



(8) 



where Qc = 7rcr Ma jcr Min is the solid angle, fi = 2.33 is the mean 
molecular weight, mu is the mass of atomic hydrogen, K Va is the 
dust opacity (or mass absorption coefficient), and B Vo (T) is the 
Planck function for dust temperature at T. We compute two dif- 
ferent column densities, one for the core A^j (with Tc and SS,) 
and the second for the total integrated flux along the line of site 
N^f 1 (with Tfuu and S™) to give an indication of the density 
of the surrounding environment. The main source of uncertainty 
here comes from the value adopted for k v . Large variations arise 
from one dust model to another, depending on the dust prop- 
erties considered : composition (with or without ice mantles), 
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Fig. 9. Distribution of the physical size of the cold clumps in pc. 
The distinction is done between the local sample (D < lkpc, 
dashed line ) and the distant sample (D > lkpc, dot-dash line). 
A power law with a - -2.3 is overlaid in dotted line over the 2 
subsets. 



structure (compact or fluffy aggregates), size... (see reviews from 
Beckwith et al. 1990; Henning et al. 1995). Dust models and ob- 
servations show that k v values can vary by a factor of 3-4 (or 
higher) from diffuse to dense and cold regions (Ossenkopf & 
Henning 1994; Kruegel & Siebenmorgen 1994; Stepnik et al. 
2003; Juvela et al. 2010). 

For this study, we have adopted the dust opacity from 
Beckwith et al. (1990) in agreement with the recommendation 
for dense clouds at intermediate densities (nn 2 ^ 10 s ) (Preibisch 
et al. 1993; Henning et al. 1995; Motte et al. 1998): 



24 



k v = O.Uv/lOOOGHz/ cm g 



2„-l 



(9) 



where we take a standard emissivity spectral index /? = 2. As vo 
is set to 857 GHz that is close to the 1000 GHz of the formula, 
the impact of variability of the spectral index ft remains small 
compared to the uncertainty of k v . For /? varying from 1 to 3, 
k Vo varies of a maximum of 15% around the value obtained with 
£ = 2. 

For the clumps for which the distance could be estimated, we 
have also determined an approximate averaged volume density 
value with : 



"h 2 =^H 2 /0"Min 



(10) 



where the third size dimension of the object is taken as equal to 
the minimum value of the clump 2D size cr Min . 

Fig. 10 shows the column densities A^ u (lower panel) and 

(upper panel) . The column density associated with the in- 
ner clump are systematically higher than the full column density, 
because it is tracing the colder and denser phase of the medium. 
We also compare the observed column density of the clump with 
Bonnor-Ebert models of cold cores (Bonnor 1956; Ebert 1955; 
Fischera & Dopita 2008) placed at 200 pc and for masses span- 
ning from 0.2 Mq (blue) to 12 Mq (red): The triangles corre- 
spond to a normal radiation field (Mathis et al. 1983) around the 
cold core, while the diamond correspond to the case of a radi- 
ation field already attenuated by external dust with Av=2. This 
modeling does not match well with the observations. One expla- 
nation is first that the Bonnor Ebert sphere modeling is only valid 
until M = 20 Mq as stressed in Montier et al. (2010), whereas 
the mass range of the C3PO catalogue is much larger as detailed 
in Sect. 4.4. Moreover it does not take into account the dilution 
inside the beam. This comparison shows also that the Planck 
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Fig. 10. Molecular column density of the clump itself iVg (up- 
per panel) and molecular column density of the total line of 
sight N^'f (lower panel) as a function of the temperature of the 
cold clump Tc- Modeling of Bonnor-Ebert spheres provides the 
temperature and column densities over-plotted in coloured sym- 
bols (triangle and diamond) for mass spanning from 0.4 (blue) 
to 12 (red) solar masses. The triangles correspond to a normal 
radiation field (Mathis et al. 1983) around the cold core, when 
the diamond correspond to the case of a radiation field already 
attenuated by external dust with Av=2. The square, cross and 
plus green symbols are respectively the very cold core in L134N 
(Pagani et al. 2004), the Pronaos core in Taurus (Stepnik et al. 
2003), and the starless cores of Herschel in Polaris Flare (Ward- 
Thompson et al. 2010). The dashed red box gives the limits of 
the domain occupied by the IRDCs of Rathborne et al. (2010). 



cold detections cannot be modeled in such a simple way and are 
probably more complex and extended objects. 

We have also over-plotted a few other reference points of 
starless cores (see caption of Fig. 10). These few objects iden- 
tified as cold cores are located in the upper distribution of the 
C3PO catalogue, in the coldest and densest part of the diagram. 
This underlines again the statistical property of the Planck ob- 
jects that have a mean column density around a few 10 21 hydro- 
gen atoms per square cm. This can be explained by the Planck 
resolution that preferentially selects quite extended objects, di- 
luting objects smaller than the 5' beam. This will be discussed in 
detail in Sect. 7. 1 . Nevertheless, we observe a few objects with 
column density greater than 10 23 , even at the Planck resolution. 
These few objects could be precursors of massive stars, or high 
mass formation regions. Moreover, the locus of the IRDCs stud- 
ied by Rathborne et al. (2010) is shown as a red dashed box in 
Fig. 10. This underlies the fact that Planck detects clumps having 
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Fig. 11. Mass spectrum for total sample (solid line), close sample 
(D < 1 kpc, dashed line) and far sample (D > 1 kpc, dot-dash 
line). A power law M~ 2 (dotted line) is overlaid for both subsets. 

the same column density but significantly colder temperatures 
than the IRDCs. 

Finally, we observe that even the densest clumps (with high 
column densities) cannot reach temperature lower than 7 K. This 
is in excellent agreement with recent observations of cold cores 
with Herschel (private communication). 

4.4. Mass Distribution 

The integrated mass over the clump is defined by: 

12 



2008), Perseus (Enoch et al. 2006), Ophiuchus (Young et al. 
2006), and Serpens (Enoch et al. 2007). They all derive mass 
function for cold cores in a range of mass 0.5 Mq < M < 20 Mq 
with slopes spanning the range a - 1.6 to 2.77 and peaking at 
around a = 2.1. An excess of high mass objects is observed 
for M > 10 4 Mq, but the heterogeneity of our sample prevents 
any further interpretation at this stage of the analysis. A better 
knowledge of the completeness of the C3PO catalogue and a 
better consistency between the distance estimates are needed. 

4.5. Luminosity 

The bolometric luminosity is defined by: 



L = 4nD 2 jSydv, 



(13) 



M 



S V0 D Z 

K V By (T) 



(11) 



where D is the distance, and S v is the integrated flux over the 
clump. The bolometric luminosity, L, is integrated over the fre- 
quency range 1 Hz < v < 1 THz, using the modeled SEDs de- 
rived from temperature and spectral index fitting (see Sect. 4.1). 
The L - M diagram is shown in Fig. 12. A large majority of the 
objects is located below the L-M line (green dot-dash) over the 
whole range of mass. The loci empirically derived by Molinari 
et al. (2008) for sources in the accretion stage (light blue) and 
in the nuclear burning stage (dark blue) are at least two order of 
magnitude above the domain covered by the C3PO clumps, in- 
dicating that accretion and nuclear burning are not dominant in 
these sources. 

The quantity L/M is very powerful in assessing the evolu- 
tionary stage of the sources, and has the advantage that it is in- 
dependent of distance: 



L ^n^Sydv 

where S Vo is the integrated flux at the frequency vo = 857GHz, D ~m~ ^ mii QcNH ' 



(14) 



is the distance, k V{) is the dust opacity (or mass absorption coef- 
ficient) as defined in Sect. 4.3, and B V(I (T) is the Planck function 
for dust temperature at T. The range of masses of the detected 
cold clumps spans from 0.3 M Q up to 2.5 x 10 4 M Q , with a me- 
dian mass of 88 M Q . The mass spectrum of the cold clumps is 
estimated by binning the mass distribution into logarithmically 
spaced bins in mass. The mass spectrum is then calculated from 



/(M) = d JL .1, 

J dM AM,- 



(12) 



where Af, is the number of clouds in bin i and AM, is the width 
of the ; th mass bin. 

As already stressed in Sect. 3.2.5, it is very difficult to char- 
acterize the completeness of the catalogue over the all-sky. The 
bias induced by the detection method inside the Galactic plane, 
due to confusion, and induced by the various methods of dis- 
tance estimate prevents any robust knowledge of the complete- 
ness of the sample. Thus the mass spectrum built here is not 
the mass spectrum of the cold core population of the entire 
Milky Way. Fig. 1 1 shows the mass spectrum of the total sample 
(solid line) and of the two subsets, D < 1 kpc (dashed line) and 
D > 1 kpc (dot-dash line), as defined in Sect. 3.2.5. A power law, 
dN/dM oc M~ a with a — 2 is overlaid (dotted line) on each mass 
function. We observe that the mass function for local objects is 
compatible with a ~ 2 over the range 30 Mq < M < 2000 Mq, 
and over the range 300 Mq < M < 10 4 M o for the distant ob- 
jects. This slope a — 2 is representative of the standard value 
a = 2.1 + 0.4 derived for MSX IRDCs with M > 100 Mq by 
Rathborne et al. (2006) . Similar mass function have been ob- 
tained on the Pipe Nebula (Alves et al. 2007) (Rathborne et al. 



Fig. 13 shows the histogram of L/M for the high-reliability 
C3PO catalogue (solid line), for the subset with distances 
(dashed line), and for all sub-samples corresponding to differ- 
ent methods of distance estimate. Three domains are defined 
following the formalism of Roy et al. (2010): Stage E (L/M < 
\Lq/Mq) corresponding to the 'Early' stage in which external 
heating is dominant ; Stage A (ILq/Mq < L/M < 30Lq/Mq) 
corresponding to the 'Accretion' -powered stage ; and the nu- 
clear burning dominant phase of the star formation (L/M > 
30Lq/Mq) . The C3PO clumps are mainly located in the L/M 
domain of Stage E for all sub-samples of the catalogue, indicat- 
ing that the nature of the objects in the catalogue is homogeneous 
with distance. Nevertheless about 15% of the C3PO sources have 
a ratio L/M > ILq/Mq and represent a candidate population 
of evolved objects in which the accretion process has already 
started and so could already contain stars. The mean tempera- 
ture of this population of clumps is around 16 K and spans from 
13 K to 18K, indicating an internal heating due to star forma- 
tion. For the other 85% of sources that fall into the early stage 
domain, it is difficult to assess the presence of stars inside the 
clumps. The lack of angular resolution prevents us from seeing 
internal sub-structures, so the presence of low-mass YSOs can- 
not be rejected (as discussed in Roy et al. (2010) for the BLAST 
population in CygX). 

Moreover we have overlaid on the L-M diagram of Fig. 12 
the curves of constant surface density £ using the theoretical for- 
mula of Krumholz (2006): 



L = 390 



M 



1 g cm- 2 100 Mq 



0.67 



Lq, 



(15) 
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Fig. 12. Bolometric luminosity as a function of mass. The L-M 
limit is over-plotted in dot-dash green line. The loci of accretion- 
powered and nuclear burning phases of Molinari et al. (2008) 
are shown in light and dark blue lines. The theoretical surface 
densities £ of Krumholz (2006) are given in dashed red lines for 
3 values: 0.01 g ■ ctrT 2 , 0.1 g ■ crrT 2 and 1 g • cm -2 . 
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Fig. 13. Histogram of the L/M ratio for the total C3PO catalogue 
(solid line) and the sub-sample of 2619 objects for which a dis- 
tance estimate has been obtained (dashed line). The distinction is 
also done between the various methods used to estimate the dis- 
tance: Molecular Complex association (green), SDSS extinction 
(light blue), 2MASS extinction (dark blue), IRDCs extinction 
(orange) and IRDCs kinematic (red). The vertical lines indicate 
the theoretical frontiers between the three stages of the star evo- 
lution: the early stage 'E', the accretion-powered stage 'A' and 
the nuclear burning dominant stage. 



for 3 values of £=0.01, 0.1 and 1 g cirT 2 (=45, 450 and 4500 
Mq pc~ 2 ). Following recent theoretical work by Krumholz & 
McKee (2008) suggesting that high-mass star form from clouds 
with S > lg ctrT 2 , it appears that only a few Planck cold clumps 
could be considered as precursors of high-mass stars, 



5. Large and medium scale distribution 

The spatial distribution of C3PO clumps is highly nonuniform; 
they seem to form arcs, groups and filaments (see Fig. 16). These 
large and small scale distribution anomalies were analysed. We 
performed an all-sky analysis (\b\ > 5°) and we show results 
for both all-sky and Tau-Aur-Per-Ori region (hereafter TAPO), 
where C3PO surface density shows remarkable excess on known 
large scale loops and shells. 



Region 




C3PO 


MC 




N G 


260 


161 ± 9 


TAPO 




0.17 


0.04 ±0.016 




6 


2.67 


2.55 ± 0.920 




No 


1833 


988 ± 25 


All-sky 




0.11 


0.06 ± 0.007 




e 


2.64 


2.54 ± 0.130 



Table 6. The number and properties of identified groups in the 
C3PO data and the Monte Carlo simulations for the TAPO and 
for the all-sky. e means the average elongation of the groups. 



5. 1 . Medium Scale Structures 

5.1.1. Groups 

We identified groups in the TAPO region using the Minimum 
Spanning Tree (MST) method of Cartwright & Whitworth 
(2004) as described in Gutermuth et al. (2009) and Beerer et al. 
(2010). A cut-off length (i.e. maximum allowed distance be- 
tween a group member core and a given subgraph) of 16 ar- 
cmins was used. It corresponds to the average distance between 
the nearest neighbours in the C3PO all-sky data. The number of 
C3PO groups, Nc in the region is 260. The fraction of groups 
with more than 3 elements, NcaINq is 17% (see Fig. 14). We 
identified groups in the all-sky data with the same method. The 
number of C3PO groups is 1833 and the value of Nga/Nq is 
11%. 

In order to assess the reliability of the statistical estimate, 
we performed a Monte-Carlo analysis consisting in 1000 ran- 
dom realizations. The same number of sources as present in the 
C3PO catalogue were randomly placed onto the sky, following 
the (l,b) marginal distributions of the C3PO positions, averaged 
in latitude and longitude bins of 5°. 

In the TAPO region in the simulated samples the average 
number of groups is 161, that is ~40% less than in the C3PO, 
see Table 6. There were at maximum 189 groups in the simulated 
samples, still 30% less then the value for C3PO. The histogram 
of the number of groups is shown in Fig. 14, where the value for 
C3PO is marked with a vertical dashed line. The all sky distri- 
bution shows similar grouping tendencies. The random samples 
show in average 988 groups, ~47% less then in the C3PO. We 
found at maximum 1083 groups in the simulated samples, still 
40% less then the value for C3PO. Our result suggests that the 
point pattern unlikely has these many groups in random distribu- 
tions. In the 1000 Monte-Carlo simulated samples the fraction of 
groups with more than 3 elements is in average 4% in the TAPO 
region and 6% in the all-sky and is never higher than 1 1 % in the 
TAPO region and 9% in the all-sky. Their histogram is shown 
in Fig. 14, where the vertical dashed line marks again the C3PO 
value. We note that the larger groups are more common in the 
C3PO than in an average random sample. 

We also investigated the variation of Nc and Nca/Nc as a 
function of the cut-off length, that we varied from 10' to 30'. We 
found that Nc and Nca/Nc increase with the cut-off length. For 
every cut-off length the average Nc and Nca/Nc in the simulated 
data are -25-45% less than in the C3PO data. 

5.1.2. Filaments 

The elongation of groups was analysed as described in Schmeja 
& Klessen (2006). We used the Cartwright & Whitworth (2004) 
definition of cluster radius: R c as the distance between the mean 
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Fig. 14. Left: histogram of the number of groups, Nc in the 
TAPO region in the Monte-Carlo simulations. Vertical dashed 
line shows the number of groups identified in the C3PO by us- 
ing the same method in the same region. Right: histogram of the 
relative number of groups with 4 or more elements Nca /Nc for 
1000 MC simulations in the TAPO region. Vertical dashed line 
shows the same value for the C3PO in this region. 
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Fig. 15. A sample group from the C3PO in the TAPO region. 
Black asterisks show the clumps, black lines indicate the MST, 
dark dashed and dotted line indicate the convex hull and the ra- 
dius of the red dashed circle is (R c ). 



position of all cluster members and the most distant sources. The 
area A of the cluster was estimated using the convex hull (the 
minimal convex set containing the set of points X in a real vec- 
tor space V) of the data points. The convex hull radius (R h ) is 
defined as the radius of a circle with an area equal to the A area 
of the convex hull of the data points. Cartwright & Whit worth 
(2004) define the elongation measure £ as follows: 



(16) 



Fig. 15 shows a sample group from the C3PO in the TAPO re- 
gion with an elongation measure of £ = 3.4. 

We calculated the elongation measure for all the 205 larger 
groups (i.e. with more than 3 members) in the C3PO all-sky data. 
We found a mean elongation of ~2.7 in the TAPO region and 
~2.6 in the all-sky, see Table 6. The mean elongation of the fila- 
ments in the Monte-Carlo simulated samples (see Sect. 5.1.1) 
does not differ from that in the C3PO in these regions, see 
Table 6. We note here the very low value of Nca in the simu- 
lated samples. We also investigated the mean elongation in the 
C3PO data for different cut-off lengths and found that the aver- 
aged elongation of the groups is insensitive to the cut-off length 
and is always around ~2.6, 



Region 






JVH_ 




IN 






TAPO 


ON 


0.912 


0.853 ± 0.023 




OFF 


0.792 


1.058 ±0.054 




IN 


0.184 


0.175 + 0.002 


All-sky 


ON 


0.322 


0.238 ± 0.005 




OFF 


0.056 


0.097 ± 0.002 



Table 7. Surface density of C3PO sources and Monte-Carlo sim- 
ulations for HI supershells in the 3 cases: IN shell, ON shell and 
OFF areas. Values for both TAPO region and all-sky are pre- 
sented. For Monte-Carlo simulations, the mean value of the dis- 
tribution is given with the 1-cr discrepancy. 



Region 




C3PO 


MC 




IN 


0.652 


1.024 ±0.033 


TAPO 


ON 


1.212 


0.906 ± 0.027 




OFF 


0.557 


0.600 ± 0.023 




IN 


0.119 


0.135 ±0.004 


All-sky 


ON 


0.193 


0.122 ±0.003 




OFF 


0.144 


0.173 ±0.002 



Table 8. Same as Table 7 for IRAS Loops 



5.2. HI shells and superloops 

A first look at the all-sky distribution of the C3PO sources 
reveals loop-like structures at medium and large scales. Here 
we quantify this impression and study the associations between 
these loops and known shells and supershells, already stressed 
to be locations of star formation. The link between the cold sub- 
structures detected by Planck and the star formation regions is 
also discussed within the framework of the triggering scenario 
(e.g. Deharveng et al. 2005; Zavagno et al. 2010a,b). 

5.2.1. HI shells 

Initially we compared the distribution of C3PO sources and the 
HI supershells from Heiles (1984). This list of shells contains 
HI shells, supershells and also shell-like objects and 'worms', 
restricted to shells crossing \b\ = 10 or located entirely outside 
\b\ = 10. These objects were originally derived from two 21 cm 
line surveys (Weaver & Williams 1973; Heiles & Habing 1974) 
exhibiting filamentary structure associated with high-velocity 
gas and radio continuum loops. Heiles (1979) observed that HI 
shells do not seem to be significantly correlated with any other 
types of object, except perhaps young stellar clusters. 

The comparison of HI shells and C3PO distribution gives the 
following results: For the TAPO region, 682 (46%) of the clumps 
are found on the shells, 635 (43%) within them and 153 (11%) 
outside the shells. For the all-sky sample the corresponding val- 
ues are as follows: 1869 (33%), 2985 (52%), 873 (15%). 

To assess the statistical reliability of these numbers, we per- 
formed an equivalent analysis on the 1000 Monte-Carlo simu- 
lations introduced in Sect. 5.1.1. The simulated surface density 
of clumps on the shells is slightly lower ( 3<x) than in the real 
sample, see Table 7. 

5.2.2. IRAS loops 

We also compared the distribution of Planck cold clumps to the 
FIR loops in the IRAS Galactic Infrared Loops of Ktiny ves et al. 
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Fig. 16. Left: Surface density grayscale map of the C3PO sources with the inner and outer boundaries of FIR loops (Kony ves et al. 
2006) overlaid. Right: Histogram of the surface density of 1000 Monte-Carlo simulated sources for ON (solid line), IN (dashed 
line) and OFF (dotted line) areas. Vertical lines show the same for the C3PO clumps, respectively. 



(2006). These IRAS loops were identified in an investigation of 
the large-scale structure of the diffuse interstellar medium (Kiss 
et al. 2004) based on 60 and 100 pm ISSA plates (IRAS Sky 
Survey Atlas Wheelock et al. 1994). They have been identified 
in the Galaxy as surfaces of high and low density ISM (Toth & 
Kiss 2007), and have been proposed as locations for star forma- 
tion (Kiss et al. 2006). These loops by definition must show an 
excess FIR intensity confined to an arc-like feature extending to 
at least 60% of a complete ellipse-shaped ring. Dust IR emis- 
sion maps by Schlegel et al. (1998) were investigated to derive 
the parameters describing the loop features (Kiss et al. 2004). 
Comparing the distribution of IRAS loops and the C3PO clumps 
gives a hint of the correlation between star formation sites and 
the cold sub-structures of the ISM. 

Fig. 16 shows the surface density map of C3PO in TAPO 
with the FIR loop boundaries overlaid. We found 810 (55%), 
312 (21%) and 348 clumps (24%) on the FIR loop shells, inside 
the loops and in the area between the FIR loops, respectively. 
Values for the all-sky study are as follows: 1928 (34%), 877 
(15%), 2922 (51%). We performed the same analysis as with 
the Heiles supershells using the Monte-Carlo simulations as de- 
scribed in Sect. 5.1.1, see Table 8. We present the histogram of 
surface density distribution of the Monte-Carlo realizations for 
the IN, ON and OFF (between loops) areas in Fig. 16b. The ver- 
tical lines give the C3PO observations. A significant (> 30%) ex- 
cess of C3PO clumps on the FIR loop shells in the TAPO region 
is found in compared to the simulations. The same behaviour is 
observed on the all-sky analysis (see Table 8). 

This analysis shows that the all-sky distribution of the C3PO 
cold clumps shows a significant correlation with known large 
scale loops, identified to be star formation regions. It is even 
more reliable in the case of IRAS loops for which this correlation 
is detected at 23 <x compared to Monte-Carlo simulations. It is 
detected at 15 cr for HI supershells on the all-sky data, but only 
at 2.6 cr in the TAPO region. 

5.3. Triggering 

In this subsection, we explore, from a qualitative point of view, 
correlations between the spatial location of Planck cold clumps 
and that of candidate sources of star formation triggering. A sig- 



nificant body of evidence (e.g. Deharveng et al. 2005; Zavagno 
et al. 2010a,b) has recently been found to support the so-called 
collect and collapse scenario, in which the swept up shell gener- 
ated by the radiation pressure associated with newly born stars is 
the preferential location for the formation of a second generation 
of stars. Alternatively, triggered star formation may occur within 
the swept up shell itself, as discussed by Elmegreen (1998). 

In this study we consider, in addition to the Planck cold 
clumps, three classes of objects, namely: the 13 most luminous 
sources recently identified by Murray & Rahman (2010) in the 
WMAP (Wilkinson Microwave Anisotropy Probe, Bennett et al. 
2003) free-free map and thought to be responsible for the bulk 
of the ionizing luminosity of the Galaxy; the 41 very massive 
star formation complexes (SFCs) found by Rahman & Murray 
(2010) to be associated with the WMAP sources; the catalogue 
of bubbles by Churchwell et al. (2006), the majority of which 
are thought to be signposts of HII regions. The ensemble of 
these three classes of objects might be representative, on dif- 
ferent angular scales, of different stages of a triggering scenario, 
with the WMAP sources being the super-giant bubbles (scale of 
multiple supermassive star formation complexes), the Rahman & 
Murray's SFCs acting as the super bubbles (scale of a supermas- 
sive star formation complex), and the Churchwell et al.'s sources 
representing the common bubbles (scale of a typical OB asso- 
ciation). Given the preliminary nature of the analysis here de- 
scribed, and due to the relatively low Planck angular resolution, 
we do not take into account smaller star formation complexes, 
such as compact and ultra-compact HII regions. 

The 13 WMAP sources are all located along the Galactic 
Plane. As for the SFCs of Rahman & Murray (2010) and 
the bubbles in the catalogue by Churchwell et al. (2006), 
these have been mostly identified using the data from the 
Spitzer GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire Benjamin et al. 2003) survey, which only covers 
the inner Galactic Plane (first and fourth quadrant). Taking into 
account the size of the WMAP sources and the coverage of the 
GLIMPSE survey, we restrict our analysis to the longitude range 
10° < 1 < 65° and 295° < 1 < 350°, \b\ < 2°. Despite the lim- 
ited area, this is where the bulk of the Galactic star formation 
takes place. At the same time, it is important to bear in mind that 
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Fig. 17. Spatial location of: Planck clumps (white solid line), Churchwell's bubbles (white dashed line), SFCs from Rahman & 
Murray (2010) (green dashed line), WMAP sources from Murray & Rahman (2010), red dashed line). Plain numbers correspond to 
entries in the C3PO catalogue, while N are names of bubbles in Churchwell's compilation. Background data are from GLIMPSE 
8/um. 
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Fig. 18. Distribution of the ECC objects (blue squares) over-plotted on the C3PO density all-sky map. 



nearby molecular clouds such as Perseus and Taurus, have not 
been included in this study. 

In the restricted coordinate region defined about, there are 
802 Planck candidate prestellar clumps. We compare the spatial 
location of these objects using, for the C3PO clumps, the size of 
the major and minor axis, and the position angles, derived from 
the fits, as described in Sect. 2.3.1 . Likewise, we make use of the 
same type of information, when provided by the authors, for the 
rest of the sources. Since at this stage we are only interested in 
a qualitative, rather than quantitative, analysis, to assess a spa- 
tial correlation we simply overlay the positions of the Planck 
sources with those of the other classes of object and count the 
number of (partial or complete) overlaps between them. In total, 
we find that 147 Planck sources (~18%) appear to be spatially 
correlated with either the WMAP sources or the SFCs (Fig. 17). 
Of these, 50 (~35%) candidate cold clumps are located at the 
edges of the super-giant or giant shells, while the remaining ob- 
jects lie within the area delineated by the shells. The relatively 
low correlation might appear surprising. However, three poten- 
tially important caveats should be noted: (1) the WMAP sources 
are visually identified in the free-free map, hence the complete- 
ness of the sample is highly uncertain; (2) likewise for the SFCs, 
given that they have been extracted with the same technique. 
In particular, the completeness will be non-uniform across the 
Galactic plane, with the more confused lines of sight being the 
most affected by incompleteness; (3) the thickness of the shells 
traced by the WMAP sources and the SFCs is not provided by 
the authors. As a consequence, many C3PO clumps which are 
currently near but not close enough to overlap with the other ob- 
jects, might actually be found to overlap when the thickness of 
the super-giant and giant bubbles is taken into account. 

We also do not find a significant spatial correlation with 
Churchwell et al.'s bubbles. However, in this case, the lack of 
correlation can be explained by the fact these bubbles and the 
C3PO sources typically have comparable angular sizes. In fact, 
often Churchwell et al.'s bubbles are smaller than the Planck can- 
didate clumps, suggesting that a triggering scenario in which the 
formation of a Planck clump is induced by the expansion of a 
GLIMPSE bubble is not very likely. 

We strongly emphasize that these results will have to be con- 
firmed by a more detailed and quantitative analysis. In particular, 
effects due to incompleteness, chance superposition and fore- 



ground/background sources along the line of sight will be inves- 
tigated using Monte Carlo simulations, and information on the 
distances of all the objects, where available. 

6. The Early Release Compact Source Catalogue 
subsample 

6. 1 . Selection Criteria 

The Early Release Compact Source Catalogue (Planck 
Collaboration 2011c) subsample, called the Early Cold Cores 
(ECC) catalogue, is a subsample of the full C3PO legacy cata- 
logue. It is obtained by applying the following two criteria: 

1. SNR > 15, 

2. r ECC < 14 K, 

where SNR is the signal-to-noise ratio of the detection, and Tecc 
is the temperature derived from aperture photometry. However, 
this temperature Tecc is derived differently from the C3PO. In 
the ECC, the photometry is done on the original Planck maps, 
i.e. maps that still retain the warm component of the diffuse dust 
emission, by placing an aperture of 5' radius on top of the de- 
tection. The background is estimated from an annulus around 
the aperture with an inner radius of 5' and an outer radius of 10'. 
Temperature is derived from a fit to all four bands used in the de- 
tection. While the annulus is efficient in subtracting the diffuse 
warm background, any warmer surrounding envelope associated 
with the cold clump remains. This estimate of the temperature 
does not take into account the real shape of the sources nor the 
subtraction of the warm envelope linked to the clump. It does, 
however, provide a simple and straightforward estimate of the 
temperature that is used only for the selection process. 

The signal-to-noise cut at SNR = 15 was chosen based on 
the Monte Carlo Quality Assessment (MCQA) results, a relia- 
bility of 90% was required, where the reliability is defined as the 
fraction of sources that have a recovered flux within 30% of that 
of the injected sources. The MCQA analysis has shown that the 
choice of SNR =15 ensures a reliability of 90% over the tem- 
perature range Tecc < 14 K in all three Planck wavelengths (cf 
Fig. 23 of Planck Collaboration 201 lv). Hence, the ECC is 90% 
reliable, where the reliability is based on flux. 

These selection criteria lead to a sub-sample of 915 objects 
over the sky and is distributed as shown in Fig. 18. Of these 
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Fig. 19. Statistical comparison between ECC and C3PO. 



objects 118 have no Simbad entry and are new detections. Thus 
at least 13% are new detections, as discussed in Sect. 2. 5. 



6.2. Comparison with C3PO catalogue 

Fig. 19 shows the histograms of various properties of the C3PO 
and ECC catalogues, normalized to the integral of each his- 
togram. The selection criteria (high SNR, low temperature) ex- 
plain the differences between the two catalogues. Because of the 
SNR selection criteria, the ECC is largely the high SNR tail of 
C3PO. The clear suppression of objects in the Galactic plane as 
shown by the latitude histogram is explained by two characteris- 
tics of the detections in this region: (1) the plane is warmer and 
hence CoCoCoDet is able to find less cold objects in the plane; 
(2) the plane is a confused region, and the SNR tends to be lower 
than for more isolated sources. The temperature shown in Fig. 19 
is not Tecc used in the selection criteria, but the derived temper- 
ature of the clump itself Tq- The ECC sources are clearly colder 
than 14 K as a rule, and even colder than the bulk of the C3PO. 

The ECC is clearly weighted towards smaller distances, 
where SNR is higher due to the flux of the cold clump, while 
more distant sources have lower SNR. Similar reasons explain 
the differences in the mass histograms; while the peak is in 
the same place, the ECC misses very low mass clumps, which 
have a poor SNR since they are faint, and misses very mas- 
sive clumps, which tend to be farther away and hence harder 
to detect. The histograms of hydrogen column density show that 
ECC is weighted towards sight-lines with higher column density, 
again due to the temperature. The longitude histograms are com- 
plex, but the peak in the direction of the anti-center is caused by 
the Taurus-Perseus-Aurigae Complex, as also seen on Fig. 18. 

Finally, the extent of the clumps, measured here by the full 
width half maximum of the fitted Gaussian profile, shows that 
the ECC sources are shifted slightly towards smaller sizes. This 
can be partly explained by the method of photometry used, but 
also by the fact that the smaller clumps tend to be colder. 

The ECC is not intended to be a complete catalogue, but a 
subset of high SNR, reliable detections of cold clumps. It is bi- 
ased towards local, cold objects which have angular sizes com- 
parable to or smaller than the Planck beam. 



7. Discussion 

7.1. Nature of the C3PO objects 

The C3PO sources have been selected by their colour properties 
as the relatively coldest Galactic sources over the sky. However 
their nature is still uncertain: are they cold clouds, clumps or 
cores? Are they really starless? The low resolution of Planck 
does not allow us to probe the inner structures of these objects to 
answer these questions. Here we will use the statistical proper- 
ties discussed in the preceeding sections to constrain the nature 
of the CP30 objects. 

The analysis of the mass-luminosity relation (see Sect. 4.5) 
has shown that 85% of the catalogue display properties of an 
early stage population, whereas the remaining 15% could be 
associated with objects that have already started the accretion- 
powered phase. Thus C3PO sources are mainly pre-stellar ob- 
jects, but not systematically pre-stellar cores. Table 9 summa- 
rizes the properties of the C3PO population in term of tem- 
perature, mass, physical size, molecular column density and 
density. Following the categorization of Williams et al. (2000) 
and Bergin & Tafalla (2007) into three populations (clouds, 
clumps and cores), the Planck detections seem mostly associ- 
ated to clumps and clouds, but not to cores that are defined to 
be smaller than 0.2 pc, with a density n > 10 4 crrr 3 and mass 
M < 5 Mq. Indeed a first study based on the prestellar core pop- 
ulation revealed by Herschel in Aquila (Konyves et al. 2010) 
and the Planck detections in the same region underlies the fact 
that Planck cold clumps are often associated with at least two 
Herschel cores. The same question is addressed in more detail 
in Paper II (Planck Collaboration 201 lr), in which a comparison 
between Planck cold clumps maps and higher resolution maps 
obtained with Herschel follow-up has been performed. 

As Planck cold clumps are tracing pre-stellar objects, a com- 
parison with IRDC properties is very instructive. A cross-match 
between C3PO sources and the MSX catalogue of 10931 IRDCs 
(Simon et al. 2006a) leads to 469 (136 inn the photometric reli- 
able catalogue) Planck Cold Clumps associated with 1007 (296 
respectively) MSX IRDCs. Notice also that less than 8% of the 
Planck clumps inside the MSX region are not directly associ- 
ated with an IRDC, or are located in the border of the spatial 
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Quantity 


nun 




max 


L c t^J 


7 


13 


17 


Nu[Hatom.cm~ ] 


10 20 


2 ■ 10 21 


2 • 10 23 


Size [pc] 


0.2 


1.2 


18 


Ellipticity 


0.4 


0.8 


1 


Mass [Mq] 


0.4 


88 


24000 


Mean density [cm -3 ] 


10 2 


2- 10 3 


10 5 



Table 9. Statistical physical properties of the Planck C3PO cat- 
alogue of cold clumps. 
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Fig. 20. Ratio of the angular extension of C3PO clumps and 
MSX IRDCs 



domain covered by the MSX survey. A comparison of the an- 
gular size of these objects is presented on Fig. 20 and clearly 
shows that the C3PO objects are more extended than the IRDCs, 
confirming their nature as clumps and not cores. Moreover the 
fact that most of the IRDCs of Simon et al. (2006a) are them- 
selves divided into cores, gives further evidence of the nature of 
the Planck detections that are only extended clumps. 

We also perform a cross-match with the Spitzer catalogue of 
IRDCs (Peretto & Fuller 2009), leading to 321 (78 on the photo- 
metric reliable catalogue) Planck Cold Clumps associated with 
1382 (356 respectively) Spitzer IRDCs. As shown on Fig. 21 
C3PO clumps could contains up to 8 MSX IRDCs and 15 Spitzer 
IRDCs. Moreover, as already stressed in Sect. 4.2, the mean ex- 
tension of the C3PO sources (about 1.4 times the Planck beam) 
is an indicator of the intrinsic scaling of r 2 for the profile of 
the extended envelope and an indicator of the physical link be- 
tween these clumps and the larger scale Galactic diatribution. All 
these observations favour the fragmentation scenario (Falgarone 
& Puget 1985; Williams et al. 2000), and place the C3PO objects 
in the middle of this process, between clouds and cores, as cold 
sub-structures. 

A detailed comparison has been performed on a smaller sam- 
ple of 14 objects for which a cross-match with MSX IRDCs has 
been obtained, and providing also kinematic distance, mass and 
density column estimates (Simon et al. 2006b). Fig. 22 shows the 
comparison of the mass, column density, density and FWHM. 
The solid line is the '1:1' relation and the dashed lines show 
the '1:10' and '1:0.1' relations. As already stressed above, the 
extension of the C3PO clumps is greater than for IRDCs. The 
C3PO clumps are 3 to 4 times more massive than the IRDCs, 
thus the mass of the cores represents only 25% to 30% of the 
total mass, as observed in the CygX cold core population (Roy 
et al. 2010). More surprising is that the column density and the 
density of the clumps are about 2 to 3 times greater than those 
of the IRDCs. This can be interpreted as a variation of the dust 



200 
150 
100 

50 




—i — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — 



MSX IRDCs 
Spitzer IRDCs 



_i . . ■ —i — ■ — ■ — l 



2 4 6 8 10 12 14 
Number of IRDCs in C3PO Clumps 

Fig. 21. Number of IRDCs found inside the Planck Cold Clumps 
when a match has been obtained with MSX catalogue (Simon 
et al. 2006a) and Spitzer catalogue (Peretto & Fuller 2009). 
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Fig. 22. Comparison of the physical properties between 14 
sources in common in the C3PO catalogue and the MSX IRDCs, 
and for which a distance estimate is available. 



opacity k v , as discussed in Sect. 4.3. We adopted the value of 
Beckwith et al. (1990) whereas it has been shown that it can 
vary by a factor of 3 or more from diffuse to dense, cold re- 
gions (Ossenkopf & Henning 1994; Kruegel & Siebenmorgen 
1994; Stepnik et al. 2003; Juvela et al. 2010). This could par- 
tially explain the discrepancy observed here. As an example, the 
Ossenkopf & Henning (1994) value of k v is about twice the es- 
timate of Beckwith et al. (1990) and would tend to lower the 
observed bias between the column densities of C3PO clumps 
and MSX IRDCs. However Simon et al. (2001) gave limita- 
tions on the column density estimates derived from CO lines 
that are known to be upper limits and to have an uncertainty 
of a few. This would then result in acceptable compatibility be- 
tween C3PO and MSX estimates inside the error bars, without 
exploring changes of the dust opacity values. 



7.2. T, /3 relation 

The variability of the dust spectral index /3 is still a contentious 
issue. Using Pronaos data Dupac et al. (2003) claimed to find 
an anti-correlation of the dust spectral index with tempera- 
ture. Since then, many investigations have been carried out 
on Archeops cold cores (Desert et al. 2008), on Boomerang 
data (Veneziani et al. 2010) and recently on Herschel data 
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Fig. 23. The relation between the temperature and spectral index 
in the catalogue and the fitted /3(T) relations. The estimates are 
based on bias corrected flux values. The relation estimated from 
IRAS and the three highest frequency HFI channels is shown 
in blue line (with lcr error in dashed line), when the red line 
is the relation obtained using only IRAS and two highest fre- 
quency HFI channels. The relation obtained on Archeops cold 
cores (Desert et al. 2008) is shown in orange. 
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Fig. 24. Monte Carlo analysis of the noise effects on the f3(T) 
relation. The black dots correspond to a synthetic sample where 
the source fluxes have been scaled to be consistent with a con- 
stant /3 value. The green points show the same sources after 
adding the estimated noise in the flux values, and the red line 
is the fit to these points. The observed relation (blue line) is sig- 
nificantly steeper. 



(Paradis et al. 2010), in which the same behaviour was found. 
Nevertheless, the issue is not straightforward because the in- 
ferred T-/3 relation suffers from a strong degeneracy between the 
two parameters. 

In the analysis described in this paper, flux values were first 
corrected for small biases estimated form Monte Carlo simula- 
tion. The derived (T,f3) values were fitted with a formula 

P = (S + coTy 1 , (17) 

(see Dupac et al. 2003). Using the IRAS 100 //m data and the 
three highest HFI frequency bands, the least squares fit gave pa- 
rameters 6 = 0.020 and co = 0.035 (see Fig. 23). When the 
analysis was repeated without the 353 GHz band, the fit give pa- 
rameters 5 = 0.053 and a> = 0.032. In particular, the presence of 
the 353 GHz band has only little effect on a), i.e. on the magni- 
tude of the temperature dependence. 



To estimate the uncertainty of the parameters, we carried out 
Monte Carlo simulations starting with the source fluxes. The un- 
certainties of the flux values were assumed to be 40% for the 
IRAS 100 jt/m and 8% for the three Planck bands. In addition, 
for each Monte Carlo realization of the catalogue, we included 
variations in the calibration. The calibration factors were gen- 
erated from normal distribution and were assumed to be 13.5% 
for the IRAS observations, 7% for the 545 GHz and 857 GHz 
and 2% for the 353 GHz Planck bands, the same scaling being 
always applied to the two highest frequency Planck bands, as 
their calibration is based on the correlation with the same FIRAS 
template (see B). With the estimated statistical uncertainties, the 
result for the analysis of the four bands is 5 — 0.020 + 0.001 and 
io = 0.035 ±0.001. 

In addition, we used Monte Carlo simulations to determine to 
what extent noise could produce an apparent correlation between 
T and /3. We started with the fluxes, temperatures, and spectral 
indices listed in the catalogue. For each source, the fluxes were 
scaled with y 2 - 2- ^ to produce a synthetic sample in which the 
spectral index is approximately constant but where the devia- 
tions from the grey body curve are still consistent with the ob- 
served sample. The resulting (T,/3) points are plotted in Fig. 24 
as black dots. In the Monte-Carlo study, noise is added to the 
fluxes, and the observed relation is no longer flat (green points 
and the red least squares line). The scatter of the points is larger 
than for the observed points (Fig. 23), suggesting that the noise 
has not been underestimated. Nevertheless, the fitted /3(T) curve 
is still flat compared to the relation observed in the catalogue 
(blue curve) and the hypothesis of a flat fi(T) relation can be 
excluded with more than 99% confidence. 

The recovered temperature dependence of the spectral index 
is much stronger than the one found in PRONAOS data (Dupac 
et al. 2003). It is also slightly steeper than the relation reported 
for Archeops cores (see Fig. 24 ; Desert et al. 2008) although the 
Archeops relations is clearly compatible with our results given 
the uncertainties of the Archeops relation. 

This kind of a f3(T) relation is expected for amorphous grains 
with disordered structure (Meny et al. 2007; Boudet et al. 2005). 
Moreover recent laboratory experiments on silica and silicates 
have shown a dependence of the spectral index on temperature 
and wavelength (Coupeaud et al. 2011a,b). As pointed out by 
Desert et al. (2008), the spectral index is expected to vary also 
with wavelength and this could have explained part of the dif- 
ference between PRONAOS at 200 - 500yum and Archeops at 
longer wavelengths, 500//m-2mm. However, with data cover- 
ing wavelengths from 100//m to 850//m, we also find a very 
steep /3{T) relation. When the spectral index is estimated with- 
out our longest wavelength band (353 GHz or 850jum) the /3(T) 
dependence becomes marginally steeper (see Fig. 23). Although 
the effect is only slight, it does suggest some flattening of the 
emission spectrum beyond 500 /mi. 

The observed colour temperature Tc and spectral index are 
not identical to the corresponding mass averaged quantities 
along the line of sight. In the presence of any observational er- 
rors, the T and j3 parameters become antic orrelated (Schwartz 
1982; Dupac et al. 2003). Shetty et al. (2009b,a) studied the 
effects of both observational noise and line of sight tempera- 
ture variations on these parameters. The colour temperature is 
biased towards the warmest regions that emit more radiation. 
Therefore, Tc overestimates the real dust temperature and the 
spectral indices are correspondingly underestimated. With a sim- 
ple model of two dust layers Shetty et al. (2009a) showed that 
Tc may rise above the real dust temperature anywhere along 
the line of sight. Similar results were obtained by Malinen et al. 
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(2010) in connection with more complex models that combined 
magneto-hydrodynamic simulations with radiative transfer mod- 
eling. Malinen et al. (2010) concluded that the spectral indices 
become more biased close to local radiation sources and this 
could affect the derived B—T relation when the sample also in- 
cludes star forming clouds. Because of all these effects, some 
caution is needed when interpreting the value and variation of 
the observed spectral index. However we see observationally 
the dust emission spectral index B increasing towards the cold 
clumps and, because temperature variations tend to decrease B 
especially towards the clump where the temperature variations 
are the largest, the actual variations in dust properties may be 
even more pronounced. 

8. Conclusion 

We have applied a dedicated source extraction algorithm, 
CoCoCoDeT (Montier et al. 2010), to the Planck data com- 
bined with the IRAS 100 pm, to build a robust catalogue of cold 
sources over the whole sky. This Cold Core Catalogue of Planck 
Objects (C3PO) is the first objectively selected all-sky catalogue 
of cold objects. It has been built using the local colour signature 
of relative colder objects embedded in a warmer background. 
We stress that this method could lead to missing sources due to 
an already cold background or spurious warm detections due to 
very hot backgrounds. The catalogue consists of 10783 objects, 
from which we have selected the 7608 sources with a complete 
set of robust physical characteristics such as fluxes, temperature, 
angular extension, ellipticity and column density. A second sub- 
sample has been constructed with 2619 sources with distance 
estimates, physical sizes, masses, and densities. In this paper we 
have performed a statistical analysis of this complete sample of 
cold objects and analysed their physical properties. 

A dedicated method has been applied to derive the photome- 
try of the clumps themselves, and to look at their local properties 
inside the warmer envelopes embedded in the Galactic environ- 
ment. The temperature of the Planck cold clumps spans from 
7K to 17 K and peaks at around 13K, in agreement with pre- 
vious studies. The advantage of the C3PO catalogue is that it 
provides a high number of very cold objects, about 600 objects 
have a temperature Tc < 10K. It has been shown that the data 
are not consistent with a constant value of the dust spectral index 
B over the whole range of temperature. Several possible scenario 
are possible, such as the effect of multiple temperature compo- 
nents folded into the measurements, and also /3(T). The question 
of the dependence of the dust spectral index B with tempera- 
ture has been discussed here and constrained, especially using 
this low temperature sample. The mean value of B is around 2. 1 
and B(T) follows a function of the temperature: B = (6 + u>T)~ l 
with 5 = 0.020 and a> = 0.035. Monte-Carlo simulations have 
demonstrated that such an anti-correlation can not be explained 
by a fitting degeneracy given our error bars. 

The mean density of these objects varies between 30crrT 3 
and 10 5 ctrT 3 with an averaged value of 2 x 10 3 cm -3 . Following 
the prescription of Williams et al. (2000), such objects are clas- 
sified as clumps, but not cores. The mass range of the catalogue 
objects varies between 0.3 and 2.5 x 10 4 Mq, the physical sizes 
span from 0.2 to 18 pc. These parameters of temperature, mass, 
density and size match well the definition of clumps given by 
Williams et al. (2000). A cross-match between C3PO sources 
and IRDCs from MSX and Spitzer has shown that each C3PO 
cold clump could contain up to 15 IRDCs, and that the C3PO 
sources are statistically more extended. Thus Planck cold detec- 
tions appear to be mainly cold clumps, intermediate structures of 



the fragmentation scenario (Falgarone & Puget 1985), between 
large clouds and very dense cold cores. Moreover, the fact that 
they are significantly extended (by a factor of almost 1 .4, com- 
pared to the PSF of the Planck-HFl instrument) indicates that 
these compact sources are linked to larger envelopes following 
a power-law profile with index in the range -2 to -1. These as- 
pects are investigated in further detail in Planck Collaboration 
(201 lr), which analyses Herschel observations of a few Planck 
sources. 

The Planck cold clump population is strongly associated 
with Galactic structures, especially the molecular component, 
and is mainly distributed within the Galactic plane. Nevertheless, 
as the detection of cold sources inside the thin Galactic plane 
is difficult due to the high confusion level, most of the de- 
tections are located in the Solar neighbourhood, within a dis- 
tance of 7kpc. By studying the correlation between the cold 
clumps and the large scale structures such as the IRAS loops 
(Konyves et al. 2007) or the WMAP triggering features (Murray 
& Rahman 2010), it has been shown in this paper that the cold 
clump population is preferably distributed on the borders of 
these large shells, where star formation is very active. Another 
interesting aspect of the C3PO is the filamentary structure of 
the clumps. This has been statistically quantified over the whole 
sky and is also apparent in the high degree of ellipticity of the 
sources. The large variety of objects covered by this all-sky cat- 
alogue gives the opportunity to perform a tentative classifica- 
tion using the physical properties, environment, and the evo- 
lution stage of these objects. In that perspective, the Herschel 
follow-up key program Galactic Cold Cores is a unique tool 
to select a representative sample of cold clumps inside this 
large C3PO catalogue, and then obtain high resolution infor- 
mation with Herschel-SPIRE/PACS. Such an analysis has al- 
ready been carried out by Juvela et al. (2010) on Herschel 
Science Demonstration Phase Data, and is continued in Planck 
Collaboration (201 lr). 

Finally, we stress that a robust sub-sample of this large C3PO 
catalogue has been delivered to the community: the Early Cold 
Core Catalogue (ECC). It provides a list of 915 objects over the 
whole sky, with similar statistics at high fluxes to the C3PO, and 
provides opportunities for ground-based or Herschel follow-up 
of cold clumps. The final legacy catalogue C3PO will benefit 
from two further Planck sky surveys. It will increase its robust- 
ness, via increased redundancy over the sky, and also increase 
the number of detections because of the reduction of the noise 
level. The polarization information will be also available in this 
final version of the catalogue. A special effort will be made to 
improve the distance estimates and so increase the volume of 
Galaxy explored, and also to improve the completeness of the 
statistical analysis. 

Appendix A: Photometry Monte-Carlo Quality 
Assessment 

As presented in Sect. 2.4, we have performed a Monte-carlo 
Analysis to assess the robustness of the photometry algorithm 
described in Sect. 2.3. A set of 10000 sources have been injected 
in Planck all-sky maps following a distribution of temperature, 
dust spectral index and flux in agreement with the observations. 
Fig. A. 1 compiles the relative errors between output and input 
quantities (fluxes, FWHM, Ellipticity, temperature and spectral 
index) for the complete set of simulated sources as a function of 
the flags raised during the photometry estimate. Three flags have 
been introduced: Aper Forced (blue), Bad Sfit 100pm (red) and 
Upper 100pm (green). 
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Fig.A.l. Comparison of the output and injected values of the Monte-Carlo simulations performed to assess the quality of the 
photometry algorithm. Black dots refer to the nominal case, when red, blue and green dots stand for various flags of the photometry 
algorithm: Bad Sfit 100pm, Aper Forced and Upper 100pm respectively. 



Appendix B: Impact of Calibration Uncertainty 

The calibration uncertainty is known to be 13.5% for the IRIS 
100 pm (Miville-Deschenes & Lagache 2005), and about 7% for 
the Planck high frequency bands at 857 GHz and 545 GHz, and 
2% at 353 GHz (Planck HFI Core Team 201 la). We investigate 
here what is the impact of this error on temperature and spec- 
tral index measurements, considering that it is not independent 
from one band to the others. The principle is to constrain the 2 
top Planck bands (545 and 857 GHz) to follow the same relative 



calibration error, when the 100 pm and the 353 GHz remain in- 
dependent. For a given temperature T and a dust spectral index 
B, SEDs are simulated using a modified black-body modeling, 
including colour correction. A 13.5% relative noise is added on 
the 100 pm fluxes, when another 7% noise realization is added 
at 545 Ghz and 857 GHz simultaneously, and a 2% noise realiza- 
tion is added on the 353 GHz fluxes. We also study the combi- 
nation of a calibration error plus an independent error band per 
band, set to 20% and 10% for 100 pm and Planck bands respec- 
tively. A third set of data is defined as our reference case and 
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Fig. B.l. Bias and lcr discrepancy of the temperature T and the dust spectral index /? estimates derived from a Monte-Carlo analysis 
dedicated to the study of the impact of the calibration error, three cases have been studied : calibration uncertainty only (red lines), 
local variance (green) and calibration plus local variance (blue). 



does not include any calibration error, but only an independent 
error per band. The T,/3 fitting algorithm described in Sect. 4. 1 
is applied on these 3 sets of SEDs, leading to an estimate of the 
bias and lcr error of the output distribution, for T ranging from 
6 K to 25 K and /? spanning from 1.75 to 2.5. 

The results are shown in Fig. B.l. The impact of the calibra- 
tion uncertainty (in red) is very low compared to the impact of 
the measurement error only (in green). In all cases, the bias on 
temperature and ft is negligible or significantly lower than 1 %. In 
both cases the bias level remains independent of T and (3. On the 
other hand the lcr error on the temperature goes from 3% to 8% 
at respectively 6 K and 25K due to calibration uncertainty only, 
and remains independent of j3. The 1-cr discrepancy on /3 is about 
10% to 7% in the temperature range and it decreases of 3% when 
P goes from 1.75 to 2.5. This error is larger and much more sen- 
sitive to /3 in the case of measurement uncertainty (green curves). 
Thus the impact of the calibration uncertainty seems very small 
compared to the impact of the measurement error of the fluxes. 
Moreover in the case of a combination of calibration and mea- 
surement uncertainties (in blue), both resulting uncertainties on 
T and /? seem to add quadratically. 

For this work, we propose to consider separately the impact 
of the calibration uncertainty, and just add the quadratic errors 
at the end of the processing. The values of the calibration uncer- 
tainties for T and f3 are provided in Annex A Table B. 1 & B.2 for 
a given set of T,/3 parameters. 
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